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I. INTRODUCTION 


A. SUPERMANEUVERABILIT Y 

In modern aerial combat, the ability of an aircraft to maneuver into the post-stall regime 
for short periods of time is crucial to its survivability. Current fighter aircraft, such as the 
Israel Aircraft Industries Lavi, SAAB Gripen or European Fighter Aircraft (EFA), employ 
a close-coupled canard to allow for continued maneuvering where conventional aircraft may 
have departed from controlled flight. This increase in maneuverability results from the 
favorable interaction of vortices over a delta wing. Double-delta wings or leading-edge 
strakes, such as those on the F/A-18 or F-16, have been used to enhance the lift in the same 
way. However, until recently, only the SAAB Viggen had been successful at using a 


canard to maintain lift at high angles of attack. 


B. CLOSE-COUPLED CANARD 

The advantages of a close-coupled canard have been known since the 1960's. It was 
found by Behrbohm [Ref. 1] that the combination of a close-coupled canard and delta- 
wing, of small aspect ratios, has significant advantages over a conventional delta-wing or 
wing/horizontal-tail configured aircraft. Both C] max and the angle of attack for CT max are 
increased by the addition of a delta-canard to a delta-wing. During the 1970's, an 
experiment was performed by Lacey [Ref. 2] to determine the correct canard geometry and 
location for maximum lift enhancement. As a result of Lacey's work, Behrbohm's 
conclusions concerning the use of a delta-canard with a delta-wing were confirmed. It was 
found that locating such a canard above, rather than coplanar with, the main wing produced 
the most favorable vortex interaction. In the 1980's canard research continued. Work by 


Er-El [Ref. 3], Stoll and Koenig [Ref. 4] and Calarese [Ref. 5] provided insight into 


canard/wing vortex interaction. However, most of their work was done at low to moderate 
angles of attack with little or no canard deflection. 

At the Naval Postgraduate School, a series of experiments has been conducted to 
compare a close-coupled canard model, designed in accordance with Lacey's work [Ref. 


2], to a wing-alone configuration. The first tests, conducted by Kersh [Ref 6], were to 


determine the forces on the model up to 50” angle of attack. During these tests, the canard 


Å o ! : _ 
was deflected between plus and minus 25 to determine which canard incidence angle 


would produced the maximum lift at a given model angle of attack. It was determined that 
the maximum lift enhancement of the canard/wing configuration over the wing-alone 
configuration occurred at 22° model angle of attack with the canard set at a positive Ú 
incidence. At this angle of attack, the first stall (loss of lift) occurred on the wing for the 
wing-alone case; the canard vortex seemed to provide a reattaching mechanism. Reference 


6 1s a complete discussion of Kersh's work. 


C. THESIS OBJECTIVE 

This thesis was the second in the series of tests conducted for canard/wing lift 
enhancement. The model was set for maximum lift enhancement conditions as described 
above. Wake surveys were then conducted using a nulling five-hole pressure probe. The 
objective was to gain quantitative data concerning the total pressure and velocity profiles at 
three crossplane locations, while the model was operating at 22° angle of attack, with and 
without the canard. By this approach the effect of the canard leading-edge vortex on the 
main wing leading-edge vortex would be further investigated. From the pressure data, 
velocity-vector plots and total-pressure contours were generated for a comparison between 
the wing-alone and canard/wing configurations. Such a comparison served to help reveal 


the enhancement mechanism at an angle of attack beyond those investigated previously. 


II. BACKGROUND 


A. VORTEX GENERATION 

The dominant characteristic of flow over a highly-swept delta wing is the generation of 
a strong leading-edge vortex as shown in Figure 1. These vortices are the result of 
separated flow at medium to high angles of attack. As opposed to the chaotic, separated 
flow associated with stall, these vortices are stable, coherent sources of high energy and 


generate areas of low Static pressure resulting in localized regions of enhanced lift at the 


leading edge. 


Primary vortex 


Secondary vortex 





Figure 1. Flowfield Over the Top of a Delta Wing 
The strongest vortices result from planforms with sharp, highly-swept (250?) leading 
edges. The sharp leading edge promotes the leading-edge separation necessary for vortex 
generation. However, the combination of large sweep and a sharp leading-edge produces a 
large planform with a low (L/D)max and a shallow lift-curve slope. The net result is poor 
range and endurance, high approach speeds and large deck space requirements. These 


characteristics are extremely undesirable traits for a carrier-based aircraft. [Ref. 2] 


The requirement for carrier suitability, therefore, dictates the design of an aircraft with 
rounded, moderately-swept leading edges. When such a planform maneuvers to high 
angles of attack the vortices produced are less coherent and of lower energy. Therefore, 
some mechanism is needed to energize or induce these vortices to remain coherent at high 


angles of attack, without the added penalty of poor cruise performance. 


B. CANARD CHARACTERISTICS 
Canards are separated into two broad categories: long-coupled or close-coupled (see 


Figure 2). 


Long-Coupled Close-Coupled 


Figure 2. Canard Categories [Ref. 2] 

Long-coupled canards are of the type used primarily as a control surface rather than as a 
lifting surface. Examples of this type are found on almost all missiles and on some aircraft, 
such as the XB-70, the Concord and the X-31 experimental aircraft. A close-coupled 
canard may provide a significant portion of the aircraft total lift in addition to being a 
control surface. The aircraft listed at the beginning of the introduction are examples of 
close-coupled designs. 

Extensive research has been performed related to the phenomenon of canard/wing 


interaction [Ref. 2,3,6,7]. These experimental results indicate two possible mechanisms by 


which the flowfield from the canard affects the flowfield around a wing. One is the 
impingement of the canard downwash on the apex of the main wing; the other is the 
favorable interference between the canard and wing leading-edge vortices. For low angles 
of attack, the canard produces a downwash field within its span and an upwash field 
outside its span. The result is a nonuniform angle of attack on the main wing. The inner 
(inboard of the canard span) and forward portions of the wing have a lower effective angle 
of attack than the outer and rear portions. The flow over the wing behind the canard tends 
to remain attached while the flow outside the canard or at the rear of the wing tends to 
separate. The probable result is an overall loss of lift on the main wing which has to be 
compensated for by the increased lift of the canard. In fact, studies such as Reference 6 
show that at low angles of attack (<10°) the lift-curve slope of a model with or without a 
canard is identical. The second mechanism occurs at higher angles of attack. As the angle 
of attack is increased, a strong leading-edge vortex is formed on the canard. As this vortex 
moves over the wing it acts to energize the wing leading-edge vortex, thereby delaying 
Vortex-Breakdown (VBD). [Ref. 3,7] 
1. Canard Versus Tail 

The first advantage of the canard over the horizontal tail design arises from a 
difference in trim requirements. A conventional tail balances an aircraft in flight by 
producing a downward lift vector. This results in an initial decrease in lift for a trim to 
lower speed. A canard, on the other hand, produces a large nose-up pitching moment 
which must be balanced by a positive elevon deflection on the main wing. The result is 
increased lift for the canard aircraft due simply to the difference in the trim requirements. 
This effect is most dramatic when comparing a pure delta-wing aircraft with a canard/delta- 
wing combination. Figure 3 shows the wing loading, approach speed and lift coefficient 


for three high-performance aircraft. The Viggen has approximately the same wing loading 


as the F-106, but the Viggen approach speed is 34 knots slower than that of the F-106 with 


A 


a 70% higher lift coefficient.[Ref. 2] 





MCDONNE LL DOUGLAS GENERAL DYNAMICS/CONVAIA SAAB 

F.4 PHANTOM (U 6) F 106 DELTA DART (U $.) VIGGEN (SWEDEN) 
WING LOADING, POUNDS PER SQUARE FOOT 64 39 40 
APPROACH SPEED, KNOTS 134 163 119 
APPROACH LIFT COEFFICIENT 1.04 0.49 0 84 


Figure 3. Characteristics of the F-4,F-106 and Viggen [Ref. 2] 
Additionally, the increased lift of a canard-configured aircraft cannot be completely 
accounted for simply by the increased lifting area of the canard. Work done by Er-El [Ref. 
3] showed an increase in normal-force coefficient of 18% at 22° angle of attack, over that 
of a corresponding wing-alone configuration. The canard itself added only 9% to the total 
lifting area. Likewise, extensive work done by Lacey [Ref. 2] showed similar results. 


Figure 4(a) shows that the stall of a wing alone and of a wing/conventional-tail combination 


o A 
occurs at about 21 , vvhereas there is no indication of stall for the canard-configured aircraft 


below 32°. Additionally, for low values of CL drag of the all three configurations was the 


same, see Figure 4(b). Therefore, during cruise flight the canard/wing configuration would 


configuration. Again, the increase in lift of the canard/wing configuration at high angles of 


attack could not be completely accounted for by the additional lifting area of the canard. 


CL, LIFT COEFFICIENT 
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Figure 4a. Lift Coefficient due to a Canard and Horizontal Tail (Ref. 2] 
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Figure 4b. Drag Polar for a Canard and Horizontal Tail [Ref. 2] 


C. VORTEX INTERACTION 

Er-E] [Ref. 3] conducted an investigation of vortex trajectory and breakdown using 
canard/wing configurations at angles of attacks between 13° and 30°. Er-El states that near 
the wing apex, the downwash of the canard forces the wing vortex down closer to the 
wing, while further downstream the vortex trajectory is more influenced by the canard 


leading-edge vortex than by its downwash. That is, near the trailing edge, the wing vortex 


is displaced upward and outward, away from the adverse pressure gradient. Er-El also 
states that a close-coupled canard causes a delay in the onset of the breakdown of the wing 
leading-edge vortex, which originates at the wing apex. This delay is possibly a result of 
the wing vortex movement away from the adverse pressure gradient, thereby making it 


more stable. 


D. PREVIOUS TESTING 

At the Naval Postgraduate School, Kersh performed force measurements on a model 
designed in accordance with Lacey's guidelines. [Ref. 6] Kersh investigated the effects of 
canard-enhanced lift at five angles of attack between 10” and 48”. Figure 5 [Ref. 6] is the 
lift curve that resulted from that study. Note that as mentioned earlier the point of maximum 
lift enhancement occurred at 22°, the point of first stall for the wing/body configuration. At 
22° there was a 34% increase in lift for the canard/wing configuration over that for the 
wing/body configuration. In addition, surface flow visualization indicated the dramatic 
effect the canard vortex had on the main wing flowfield. However, while the flow 
visualization gave outstanding insight into the nature of the flow on the surface of the wing, 
it did not say much about what occurred off the wing surface. That flowfield study was the 


major objective of this study. 
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Figure 5. Lift Coefficient of Wing/body and Maximum Lift of Canard/Wing 
[Ref. 6] 
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HI. EXPERIMENT AND PROCEDURES 


A. PURPOSE 


This study was conducted to Investigate the effects of vortex interference on the 
flowfield around the main wing of a close-coupled canard and vving combination. A nulling 
five-hole pressure probe was used to make flowfield pressure measurements in the three 
crossplanar grid locations shown in Figure 6. These pressure measurements were then 
used to determine local flow velocities, pressures, and pitch and yaw angles for use in 


mapping the position and interaction of canard and wing vortices. 


Platform 


Freestream Direction 


TIO DE 
: 35 xO" (3) 
4"x5" (2) 
(1) 


Measuring grids 





Figure 6. Model Overview 


B. APPARATUS 

The primary equipment used in this study vvere a vvind tunnel, a close-coupled canard 
model, a three-dimensional traversing mechanism, a rotary pressure transducer, a data 
acquisition system, and a nulling five-hole pressure probe. 

1. Wind Tunnel 

The tunnel was of the close-circuit, single-return type. It measured 64 feet in length 

and between 21.5 and 25.5 feet in width with a test section cross-sectional area of 8.75 
square feet. The tunnel was powered by a 100hp electric motor coupled to a three-bladed 


variable pitch fan via a four-speed transmission. Figure 7 is a schematic of the tunnel. 
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Figure 7. Naval Postgraduate School Wind Tunnel [Ref. 8] 
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Tunnel turbulence was minimized by the use of stator blades, turning vanes and 
turbulence screens. Eight stator blades were used to remove swirl imparted to the flow by 
the fan blades. Turning vanes were located at each corner to turn the flow while 
maintaining its kinetic energy and turbulence screens were placed in the settling chamber. 

The test section measured 45 inches by 28.5 inches with a cross-sectional area 
about one tenth that of the settling chamber. The test section was rectangular with slightly 
divergent walls to account for boundary layer growth. Breather slots at the far end of the 
test section were used to keep the test section at approximately atmospheric pressure. Due 
to the configuration of the pressure measuring equipment, ambient atmospheric pressure 
was used as the reference static pressure rather than test section static pressure. [Ref. 8] 

Tunnel velocity was set by reference to a water manometer that measured the 
difference in static pressure between the settling chamber and the test section. The static 
pressures in the settling chamber and at the test section entrance were determined by 
reference to four pressure ports, one on each wall. The pressures from these ports were 
averaged by a common manifold prior to the manometer. The manometer measured the 
pressure difference in centimeters of water and the test section velocity was then determined 


by equation (1). 
V = {(2*2.0475*PcmH20)/p*K)} 1/2 (1) 


where: 
V Test Section Velocity (ft/sec) 
2.0475 Converts From Centimeters of Water to Lb¢ / ft? 
PemH20 Manometer Reading (cm of water) 
p Density of Air (slugs/ft3) 


K Tunnel Calibration Factor (PemH20/q) 
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The tunnel calibration factor, the ratio of the static-pressure ratio to the test section 

dynamic pressure, was determined to be 0.925. 
2. Model Design 

The model used in this study was originally used by Kersh [Ref. 6] and was 
designed in accordance with the earlier work of Lacey [Ref. 2]. During Lacey's studies, a 
canard-area-to-wing-area ratio of 0.20 as referenced to the model centerline of the fuselage 
was used. Because it was desired to accommodate an internal mechanism for a movable 
canard, the fuselage of this model was larger than that of Lacey's model. Therefore, this 
model was designed with an exposed canard-area-to-wing-area (wing referenced to the 
centerline of the fuselage) ratio of 0.20. If the area ratio had not been determined this way, 
a canard exposed area much smaller than Lacey's would have resulted. If the area ratio of 
Lacey's model had been determined by this same method, a ratio of 0.13 would have 
resulted. Consequently, an exposed canard-area-to-wing-area ratio of 0.20 compared 
favorably with that of Lacey's model. 

Lacey's work also showed that the relative position of the canard to the wing was 
crucial for constructive canard/wing interference. According to Reference 2, the 
longitudinal position of the canard (X), non-dimensionalized by the mean aerodynamic 
chord (MAC) of the wing referenced to the fuselage centerline (X/Cmac), should not be 
greater than X/Cmac=1.5. The longitudinal distance was measured from the 40% chord 
location of the exposed canard root to the quarter chord of the wing MAC. Also, there 
should be no overlap of the main wing and canard. As a result, the canard was positioned 
at X/Cmac=1.2. The canard vertical position was Z/Cmac=0.2; which was chosen such 
that it would match Lacey's model. 


The wing and canard were both an NACA 64A008 airfoil section, which was the 


same section used by Lacey. The canard had a leading-edge sweep of 60”, a straight-taper 


of 0.1 and an aspect ratio of 2.0. The wing had a leading-edge sweep of 50°, a straight- 


taper of 0.15 and an aspect ratio of 3.0. These planforms vvere chosen to ensure that strong 
leading-edge vortices would be generated. Unlike the sharp leading-edge model used in Er- 
El's studies [Ref. 3], both the canard and wing had rounded leading edges. No attempt was 
made to trip the boundary layer. Planform geometry was derived from equations (2), (3) 


and (4). [Ref. 6] 


AR z 2(b/C,(14A)) (2) 

AR - b2/S (3) 

MAC z 2/3(Cr + Ct - [CrC//(Cr + Ct)]) (4) 
Where: 


AR Aspect Ratio 

b Span 

Cr Length of Root Chord 

Ct Length of Tip Chord 

À Taper Ratio (C,/Cr) 

S Area of VVing or Canard 

MAC Wing Mean Aerodynamic Chord 


Figure 8 is a schematic of the model and contains a listing of its geometric data. 
Reference 6 contains a thorough description of the model's design. Because the five-hole 
pressure probe was not long enough to reach the model mounted on the tunnel floor, a 
platform was constructed to elevate the model (see Figure 9). The effect on the flowfield 
should be minimal for the half-model, as the raised platform acted in a similar manner to the 


reflection plane. 
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MAC=7.55" 





Airfoil Section (NACA) 64A008 64A008 
Area (semi-span), in” 

Projected Area 97.9 50.7 

Exposed Area 59.6 1953 
Exposed Semi-Span, inches 9.1 4.4 
Centerline Semi-Span, 123] 7.4 
inches 
Chord, inches 

Exposed Root 11.0 8.0 

Centerline Root 14.0 b 

Tip 21 0.8 
Aspect Ratio 3 2 
Taper Ratio 0.15 0.1 
Sweep Angle in Degrees 

Leading Edge 50 60 

Trailing Edge 10.6 ə. 
Incidence Angle in Degrees 0 - 
Dihedral 0 0 
Twist 0 — 


Figure 8. Model Geometric Data 
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Figure 9(b). Front Vievv of the Model Mounted in the Tunnel 
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3. Velmex 8300 Three-Dimensional Traverser System 

The Velmex system, shovvn in Figure 10, vvas used to accurately position the five- 
hole pressure probe in the test section. The system consisted of the traverser mechanism, 
Figure 10(b), and the motor controller, Figure 10(a). The traverser mechanism vvas 
composed of leadscrews, slides and motor/jackscrew assemblies. Each motor was a 200- 
step-per-revolution, 10-amp stepping motor with a maximum velocity of 3000 steps per 
second. The motor step size was 1/200 of a revolution which equated to 0.000125 inch. 
The three motor assemblies received their commands from the controller which was capable 
of interpreting signals from either a parent computer or manual inputs at the controller. 
During this study the controller received its commands from a parent PC/AT.[Ref. 9] 

4. Rotary Pressure Transducer and Data Acquisition System 

A 48-port rotary pressure transducer (commonly referred to as a Scanivalve®) read 
each port of the pressure probe twice per measurement point. The data acquisition system, 
shown in Figure 11, consisted of the required hardware and software to allow a PC/AT to 
acquire the pressure data. The Scanivalve® put out a 7-bit binary coded decimal signal that 
corresponded to the port currently being monitored. This arrangement allowed remote, 
electronic monitoring and control of the Scanivalve® assembly. [Ref. 10] After the 
pressure signal, in the form of an analog voltage, left the Scanivalve® it was passed to the 
Relay Multiplexer, then to the Digital Multimeter (DMM). The DMM converted the analog 
voltage signals to a din form which were then sent to the computer for further 


processing. 
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Figure 10(a). Velmex 8300 Traverser - Motor Controller 
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Figure 10(b). Velmex 8300 Traverser Assembly 
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Figure 11. Data Acquisition Hardware [Ref 10] 
5. Five-Hole Pressure Probe 
The three-dimensional five-hole pressure probe, shown in Figure 12, was made of 
corrosion-resistant, non-magnetic, stainless steel. It was 0.125 inches in diameter and 24 


inches in total length. At the tip of the probe were five measuring holes. The center hole of 


the probe (P4) measured total pressure, the two lateral holes (P2,Pz) were used to measure 


yaw angle and static pressure, and the top and bottom holes (Pa,Ps) were used for pitch 


information. [Ref. 11] 
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Figure 12. The Five-Hole Pressure Probe [Ref. 11] 
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C. EXPERIMENT SOFTWARE 
The software used in this study was developed by Lung during his work with the five- 
hole pressure probe [Ref. 10]. Some modifications were made to accommodate this work, 
but the flow of data remained essentially the same. Figure 13 is a schematic flowchart of 
the various programs and data files. The following sections discuss the individual 
programs. 
1. PPROBE 

PPROBE (Appendix A), written in BASICA by Kindelspire [Ref. 12], and 
modified by Lung [Ref. 10], was used to control the data acquisition process. Figure 14 is 
the program flow chart for PPROBE. The two primary responsibilities of PPROBE were 
the control of the traverser and the recording of the data. 

First, PPROBE was used to manually position the probe to the proper starting point 
in the tunnel, then the dimensions of the grid to be measured and the desired step size were 
entered. PPROBE calculated the total number of points to be measured, created the required 
number of data files and initialized the data acquisition equipment. The operator then 
physically nulled the probe and entered the observed yaw angle into PPROBE. PPROBE 
rotated the pressure transducer to the appropriate port and began measuring the analog 
voltage which corresponded to the pressure at that port. Finally, the data were displayed to 
the operator and he was given a chance to remeasure the data point. If the data were 
acceptable PPROBE moved the traverser to the next measuring point. At the completion of 


each column in the grid, the data were placed into a file to be used in the program PVA. 
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Figure 13. Data Flow Path 
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Figure 14. PPROBE Flow Chart 
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2. CALP 
CALP (Appendix B), was written in BASICA and used to calibrate the pressure 
transducer prior to each data collection period. The calibration manometer is shown in 
Figure 15. Six pressures were recorded using the calibration manometer and the pressure 
transducer in the Scanivalve® to provide a known calibration curve for the experiment. A 


straight line was fit to the data, whose slope and intercept were input to the PVA program. 
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iF OTHER 
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Figure 15. Calibration Manometer 
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3. PVA (Pressure Velocity Angle) 

PVA (Appendix C) was a modified version of CONVERT, a Fortran program 
which combined the output of PPROBE and CALP to calculate the flowfield properties 
[Ref 10]. PVA used the pressure data, in the form of analog voltages, and yaw data from 
PPROBE, and combined them with the curve fit derived from the CALP output to 
determine the desired flowfield characteristics. See Figure 16 for the PVA flow chart. 

Specifically, PVA determined: 

-freestream velocity (V) 

-crossflow components of freestream velocity (Vx and V y) 
-total pressure and total pressure coefficients 

-Static pressure and static pressure coefficients 

-flow yaw and pitch angles 

The pressure coefficients were determined as follows: 

CPT = {PT1-Pamb-Q1}/Q1 (5) 

Cps = {Ps1- Pamb}/Q1 (6) 

Where: 

PT] = Local Dynamic Pressure 

Ps] = Local Static Pressure 

Pamb = Room Ambient Pressure 

Qı =Freestream (Tunnel) Dynamic Pressure 
Due to the configuration of the pressure-measuring equipment, room ambient pressure was 
used as the reference static pressure in the pressure coefficient calculations. The output of 
PVA was put in an ASCII file for use by commercially available graphic software for 


plotting the desired output. 
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PVA Flow Chart (1) 


Input 
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Figure 16(a). PVA Program Flow Chart 
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From PVA (1) 


Calculate Q, Q1 
Q =local dynamic pressure 
Q1=freestream dynamic pressure 
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Figure 16(b). PVA Program Flow Chart 
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D. EXPERIMENTAL CONDITIONS AND PROCEDURES 


To correlate these results with the results of Kersh, the experimental conditions were 


kept as constant as possible. A reference measured pressure difference of 17 cm of H20 
was Set, resulting in a test-section dynamic pressure of 37.63 Ibg/ft? and an average 
velocity of 177 ft/s. The Reynolds number was Re=9.5x10° based on the main wing mean 
aerodynamic chord. 

To account for possible tunnel crossflow, probe installation effects and flowfield 
disturbances caused by the model platform, an initial run was made without the model 
installed, to established baseline yaw and pitch values. This measurement was again run at 
the end of testing to ensure the baseline values had not changed appreciably. 

Measurements were made at three locations, numbered as shown in Figure 6. Grids 1 
and 2 were mapped at 1/4" intervals, while grid 3 was mapped at 1/2" intervals. Note that 
the measuring plane of the grids were perpendicular to the freestream direction, rather than 
to the wing chord. 

Grid 1 was 4" by 5" and was set just behind the canard trailing edge in order to capture 
the canard leading-edge vortex. Grid number 2 was 5.5" by 6" and was placed at the mid- 
point of the wing to capture the effect of the canard leading-edge vortex on the flowfield at 
that point. Finally, since flow visualization, completed earlier, had shown the flowfield 
beginning to lose its coherency near the trailing edge of the wing, grid number 3 was 
placed just behind the trailing edge of the wing. Grid 3 measured 7.5" by 6"; because its 
size would have necessitated an excessive amount of tunnel time to complete at 1/4" 


intervals, it was mapped at 1/2" intervals. 
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IV. RESULTS 


The following sections discuss the results of this study. In each section there will be 
three plots. The first will be a plot of flowfield velocity vectors, then a plot of the velocity 
streamlines, followed finally by a plot of total-pressure-coefficient contours. Flow 
visualization results are also included to correlate with the measured data. All the plots are 
scaled in inches, with "X" plotted horizontally, transversely to the tunnel and "Y" plotted 
vertically to the tunnel. Locations on the plots will be identified by listing the coordinates in 
the standard (X,Y) convention. Also, the plotted data are shown with cutaway views of the 
model geometry to establish the relative position of the measuring plane to the model. At 
grid 1 the actual measuring plane was just beyond the canard trailing edge, but the canard 
cross-section at the trailing edge will be shown for reference. Likewise, at station three the 
grid was located after the main wing trailing edge, but the wing cross section 3 is shown to 


indicate the relative position of the vortices and main wing (see Figure 17). 





Figure 17. Model Side view 
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Because of the dynamic nature of flowfields surrounding aerodynamic surfaces at high 
angles of attack, there were regions in which the capabilities of the probe were exceeded. In 
particular, the probe was unable to make accurate measurements downstream of stalled 
aerodynamic surfaces or near the core of a vortex. Precisely, in these regions the probe 
could not be nulled - that is, P2 and P3 could not be made equal. The main effect of this 
probe limitation was that the exact positions of the vortex cores could not be determined. 
During data reduction, to account for the regions of the flowfield where the probe could not 
be nulled, a value of zero velocity was assigned and a limiting value of total pressure drop 
was established based on observed trends in the area of the flowfield where accurate data 
were available. Therefore, the velocity vector plots and velocity streamline plots have areas 
where there are no data present; these areas are the sections of the flowfield where the 
probe capabilities were exceeded. Likewise, on the pressure contour plots an area of 
uniform total pressure loss indicates the section in which the data are inaccurate. The isobar 
along which the pressure data become unreliable will be delineated in the section which 
refers to that plot. 

An additional problem of the probe's lack of capabilities arises from the fact that both 
vortical flow and the reversed flow from a stalled aerodynamic surface appear as blank 
areas in the data or as areas of total pressure loss. Since there is insufficient data to indicate 
a difference between the two flow states, flow visualization results have been included to 
discern between regions of separated/reversed flow and those of separated coherent vortex 


flow. 


A. WING-ALONE 
With the canard off, the flow over the main wing was characterized by large regions of 
reversed flow, indicative of a stalled wing. In general, as can be seen from Figure 18, there 


was a large amount of oil over the majority of the wing indicating wing stall. Particularly, 
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the oil on the inboard portion of the wing was mostly stagnate or flowing spanwise; this 
response indicated that this section of the wing had stalled and thereby created an area of 
large total pressure loss. Along the leading edge of the wing there appeared to be a weak, 
poorly organized, separated vortex, which apparently lost its coherency over the outboard 
section. This separated vortex was indicated in Figure 18 by the heavy line of oil which 
started at the wing apex and then seemingly burst at the wing mid-point. Therefore, it 
would be expected that at each measuring grid there would be evident a general tendency 
for the flow to move spanwise, in a root-to-tip fashion, with a large area of total pressure 
loss. In this case, the loss in total pressure would be due to chaotic, disorganized, 


separated reversed flow rather than coherent, vortical flow. 





Figure 18. Wing Alone Set at 27 
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1. Mid-Point of the Main VVing (Grid.2) 

Figure 19(a) is a plot of the "X" and "Y" components of the local velocities at gnd 
2. The blank area inside the heavy line indicates the area where the data were unreliable due 
to flow separation. The crossflow component of freestream velocity increased as the probe 
was traversed from the wing tip to the root. The crossplane velocity streamline plot of 
Figure 19(b) shows a vortical pattern which compares favorably with the flow visualization 
of Figure 18. Crossplane velocities are vertical (in the model reference sense) near the 
leading edge or tip, but horizontal at 26% mean aerodynamic chord (MAC) above the wing. 
Figure 20 is a contour plot of total-pressure-coefficient loss. This plot also shows an area 
of total pressure loss increasing in size as the probe was moved inboard toward the wing 
root. At the outboard section of the wing it is possible that the cause of this loss was the 
disintegrating vortex that appeared to form there; however, at the inboard section, the loss 
was certainly caused by reversed flow. Pressures inside the -1.24 isobar are unreliable. 
The separated region extends up to 29% MAC. 

2. Trailing Edge of the Main Wing (Grid 3) 

At grid 3, which was located just past the trailing edge of the wing, Figures 21(a) 
and (b) clearly indicated a strong spanwise flow pattern from the wing root to the tip near 
the fuselage. This flow pattern is believed to be caused by entrainment due to the separated 
reversed flow commonly associated with wing stall, rather than by a coherent vortex. The 
separation now extends about 57% of the MAC above the wing plane at the trailing edge. 
The total-pressure-coefficient gradients in Figure 22 show that the region of total pressure 
loss, due to the stalled reversed flow, had grow substantially as the flow moved 
downstream from grid 2. Pressures inside the -1.20 isobar are considered to be unreliable. 
Though the region of uncertain data is large, the results seem to indicate the extent to which 
the flow separation affects the wing flowfield, just as the oil flow visualization seemed to 


indicate the same on the wing surface. 
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B. VVING/CANARD COMBINATION 

The addition of the canard to the model caused large-scale reattachment of the flovv over 
the wing. Figure 23 shows a well-defined secondary separation line along the leading 
edges of both the canard and wing. These separation lines are characteristic of strong, 


highly coherent leading-edge vortices. Note in Figure 23 that despite the high angle of 


attack of the canard (29°), the flow appeared to remain completely attached (not separated 
and reversed) and that the large leading-edge sweep (60°) produced a strong vortex even 
though the leading edge was rounded. In reference to Figure 23, the wing can be seen to be 
fairly clean of oil with the exception of a wedge which started just aft of the apex and grew 
as it proceeded downstream parallel to the fuselage. Also, there is a line of oil along the 
trailing edge of the wing which indicated that the flow had begun to separate there, 
particularly at the point were the trailing-edge line joined the line which had started near the 
apex and moved parallel to the fuselage. Additionally, looking at the top surface of the 
fuselage where the canard and wing join the body, one can see there is a section of the 
fuselage free of oil. It is thought that this result was due to a secondary vortex forming near 
the surface of the fuselage. The exact origin of this secondary fuselage vortex, however, is 
not clear. It could have been the result of either the forebody or of the leading-edge vortices 
formed on the canard and wing, or possibly some combination of the two effects. 

The discussions that follow will be in the same format used in the preceeding wing- 
alone discussion. The path of the canard and wing vortices will be analyzed and the effects 


of the one on the other noted. 
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Figure 23. Model set at 22 Canard set at 7 Incidence 


1. Trailing Edge of the Canard (Grid 1) 

The flow at this grid, which was located just aft of the canard trailing edge, was 
dominated by a strong vortex generated at the canard leading edge. Figure 24(a) shows that 
for a large portion of the measuring grid, the probe was unable to accurately record flow 
field pressures. However, Figures 24(a) and (b) do indicate the presence of very strong 
vortical flow in a clockwise direction. Also, Figures 24(a) and (b) indicate the presence of a 
secondary vortex forming toward the fuselage surface near the canard/body juncture. The 
secondary vortex was rotating in a counter-clockwise direction, which indicated that it had 


probably formed as a result of the canard leading-edge vortex in much the same manner as 
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a secondary vortex forms near a vving leading-edge vortex. The mechanism for such a 
vortex would be the same as that depicted in Figure 1, but would be a result of an 
interaction of the wing and fuselage rather than the wing and itself. Such a phenomenon 
had not been noted in previous testing because the models used did not have a substantial 
fuselage. If this secondary vortex had resulted from the model forebody, it would have 
been expected to be rotating clockwise. 

Figure 25 shows the steep pressure gradients which resulted from this vortex. The 
data inside the -1.26 isobar are unreliable, but clearly, the trend was for a tight, coherent 
core vortex standing approximately 1.5 to 2.5 inches off the canard surface at the trailing 
edge. The presence of the secondary vortex near the canard root is better indicated in the 
plot of total-pressure-coefficient-loss contours. Notice that a freestream value of total 
pressure extends all the way to the canard surface, which agrees with the flow visualization 


of Figure 23. 
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2. Mid-Point of the Main VVing (Grid 2) 

A comparison of Figures 26 and 19 shows a dramatic difference in the flowfield 
due to the addition of the canard to the model. Figure 26(a) shows two regions of 
unobtainable data, one at approximately (5,2) and the other very near the wing leading edge 
at about (0.5,3.5). The area located near (5,2) was the canard leading-edge vortex which 
had moved downstream from grid 1. The area near (0.5,3.5) could have been the result of 
two possible phenomena. One is that the large region of chaotic flow seen in Figure 19 had 
been forced into a small "pocket" of separated reversed flow. The second is that the wing 
leading-edge vortex had reformed due to the presence of the canard vortex. The flow 
visualization in Figure 27 indicates the presence of a wing leading-edge vortex. 
Furthermore, the pattern formed by the streamlines in Figure 26(b) indicates that rather than 
the flow having merely been forced into a "separation-pocket", the leading-edge vortex was 
actually energized by the canard vortex and had reformed. The flow visualization suggests 
that a strong vortex has formed along the wing surface. 

The energizing of the wing leading-edge vortex in this manner and the reattached 
flow over the inboard wing section would lead to a level of enhanced lift which would be 
greater than that obtained by simply adding the lifting surface of the canard. In fact, as was 
mentioned earlier, this was found to be exactly the case in previous studies [Ref. 2,3]. The 
mechanism is clearly demonstrated here. Also in Figure 26 can be seen a significant 
crossflow component between the canard vortex and the wing leading-edge vortex. In the 
vicinity of (2,1) the crossflow is as much as 65% of the total freestream velocity and 
approaches the wing perpendicular to its surface. The formation of oil parallel to the 
fuselage in Figure 27 could have been due, at least in part, to the stagnation of this large 
crossflow component on the wing surface. 

The total-pressure-coefficient contours of Figure 28 show the presence of a third 


vortex not completely captured by the velocity-vector or velocity-streamline plots. Centered 
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at (2.1,0.6), this vortex is believed to be the secondary vortex that had formed on the 
fuselage surface at grid 1. Most importantly, the contours of Figure 28 show that 
freestream total pressure extended to the wing surface over the inboard section. Clear 
evidence is provided that the flow had reattached over the inboard portion of the wing, 


enhancing significantly the lifting ability of the previously stalled wing. 
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Figure 27. VVing at 22° with Canard Deflected 7° 


49 





PAEUED Ҹилд ‘SULM JO jutod-prA 


"SINOJUOJ JUJVIJJJOL) ә1П559ЈД 19301, *82 


JINŻIH 





50 


3. Trailing Edge of the Main VVing (Grid 3) 

VVhile the velocity vectors and streamlines of Figure 29 indicate a similar pattern to 
that of grid 2, the flow, in general, was less coherent and slightly more diffuse than at grid 
2: 

Figure 29(a) shows that there were two areas devoid of obtainable data. One of 
these areas was the result of the canard vortex which had moved downstream to a point 
near the (5.5,1.5) position and the other was caused by the wing leading-edge vortex 
which had reformed as a result of the presence of the canard. The streamlines in Figure 
29(b) show that the flow at these locations was strongly vortical. In addition, the flow 
visualization of Figure 27 shows the presence of a secondary separation line consistent 
with the formation of a leading-edge vortex. Figure 29 also indicated that there still existed 
a strong crossflow component moving from the canard vortex to the wing vortex. The 
magnitude of this crossflow reached 66% of the freestream as at grid 2; however, the 
freestream velocity on the wing surface of grid 3, near (0,2.5), had slowed to 74% of that 
at the same position on grid 2. This was substantiated by the growth in the width of the oil 
line parallel to the fuselage. A comparison of Figures 18 and 21 with 29 shows that the 
energy introduced to the flow by the canard vortex moved the burst point of the wing 
leading-edge vortex past the trailing edge of the wing. 

The pressure contours in Figure 30 show the presence of the two vortices already 
mentioned and a third area of pressure loss. This third area was most likely a combination 
of the remnants of the secondary fuselage vortex and the stagnation of the spanwise 
crossflow component on the wing. Note also that even though the canard vortex had 
traveled the length of the model, the pressure gradient around it was approximately equal to 
that around the wing leading-edge vortex. This further indicates that it was the energy in the 
canard vortex which allowed or entrained the wing vortex to reform. The -1.21 isobar 


marks the limit of the probe capabilities. 
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C. CONCLUSION 

The addition of the canard had a dramatic effect on the wing flowfield at an angle of 
attack of 22°. Without the canard, the flow across the wing was characterized by large 
areas of chaotic, incoherent flow with a resultant loss in lift. The addition of the canard 
vortex established a large crossflow component moving toward the wing surface which 
caused the flow over the inboard section of the wing to reattach. 

In addition to reattachment, the crossflow provided energy to the wing vortex which 
delayed its breakdown. Therefore, with the canard in this position relative to the wing, the 
the wing leading-edge vortex was enhanced and stabilized. This crossflow-induced 
breakdown delay substantiates the conclusions made by Lacey, in Reference 2. 

The energy contained in the canard vortex is not lost. Consider the drag force as a 
measure of the energy required to propel two models with the same wing planform. The 
first model employs a canard/wing configuration and the second model a 
wing/conventional-tail configuration. If the flow over the wing of the second model 
remains attached (this study shows it would not), both models would have approximately 
the same drag. The difference is that the energy (in the form of drag) placed in the tail 
vortex of the second model would be lost to the freestream. In the canard/wing 
configuration, some of the energy placed in the canard vortex was used to benefit the flow 
over the wing thereby, making the canard configuration more efficient. 

Figure 31 1s provided as a schematic of the approximate path of the vortices as they 


moved downstream through the grids. 
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Freestream Direction 





Figure 31. Vortex Path 
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APPENDIX A 


1 DEF SEG : CLEAR , &HFEOO: GOTO 4’Begin PCIB Program Shell 
2 GOTO 1000)” User progran 

3 GOTO 900 * Error handling 

4 I = &HFEOO’ Copyright Hewlett-Packard 1984,1985 
5 PCIB.DIR$ = ENVIRONS ("PCIB") 

6 IS = PCIB. DIRS + "\PCIBILC. BLD" 

7 BLOAD I$, I 

8 CALL I(PCIB.DIRS, 1%, (02): PCIB.SEG = 15 

9 TE Js = 0 THEN GOTO 13 

10 PRINT "Unable to load."; 

11 PRINT " (Error or". Uši mi 

12 END 

1321 


14 DEF SEG = PCIB.SEG: O.S = 5: C.S = 10: I.V = 15 

15 I.C = 20: L.P = 25: LD.FILE = 30 

16 GET.MEM = 35: L.S = 40: PANELS = 45: DEF.ERR = 50 

17 PCIB.ERRS = STRINGS(64, 32): PCIB.NAMES - STRINGS(16, 32) 
18 CALL DEF.ERR(PCIB.ERR, PCIB.ERR$, PCIB.NAME$, PCIB.GLBERR): 
PCIB.BASERR = 255 

19 ON ERROR GOTO 3 

20 J = =1 

21 I$ = PCIB.DIR$ + "\PCIB.SYN" 

22T OSS) 

23 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

24 Т = О 

25 CAME  I.V(I, READ. REGISTER, READ.SELFID, DEFINE, 
INITIALIZE.SYSTEM) 

26 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

27 CALL I.V(I, ENABLE.SYSTEM, DISABLE.SYSTEM, INITIALIZE, 
POWER. ON) 

28 IF PCIB.ERR «» 0 THEN ERROR PCIB.BASERR 

29 CALL I.V(I,MEASURE, OUTPUT, START, HALT) 

30 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

31 CALL I.V(I, ENABLE.INT.TRIGGER, DISABLE.INT.TRIGGER, 
ENABLE.OUTPUT, DISABLE.OUTPUT) 

32 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

33 CALL I.V(I,  CHECK.DONE, GET.STATUS,  SET.FUNCTION, 
SET.RANGE) 

34 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

35 CALL I.V(I, SET.MODE, WRITE.CAL, READ.CAL, STORE.CAL) 

36 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

37 CALL I.V(I, DELAY, SAVE.SYSTEM, J, J) 

38 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

Son jil 

40 CALL  I.V(I,  SET.GATETIME,  SET.SAMPLES,  SET.SLOPE, 
SET. SOURCE) 
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IF PCIB.ERR 
ee I.C(I, 
IF PCIB.ERR 
2 EL I.cÇI, 
IF PCIB.ERR 
SMEL I.C(I, 
IF PCIB.ERR 
SAME I.C(I, 
IF PCIB.ERR 
BALL I.C(I, 
IF PCIB.ERR 
I - 2 

m-ə 

7 I.V(I, 
IF PCIB.ERR 
SEM I.C(I, 
IF PCIB.ERR 
RE I.C(I, 
IF PCIB.ERR 
EXE I.C(I, 
IF PCIB.ERR 
EXER T.C(I, 
IF PCIB.ERR 
I = 4 

EXER I.V(I, 


«» 0 THEN ERROR PCIB.BASERR 
FREQUENCY, AUTO.FREQ, PERIOD, AUTO.PER) 
«» 0 THEN ERROR PCIB.BASERR 

INTERVAL, RATIO, TOTALIZE, R100MILLI) 
«» 0 THEN ERROR PCIB.BASERR 

R1, R10, R100, RIKILO) 

«» 0 THEN ERROR PCIB.BASERR 

RIOMEGA, RIOOMEGA, CHAN.A, CHAN.B) 

«» 0 THEN ERROR PCIB.BASERR 

POSITIVE, NEGATIVE, COMN, SEPARATE) 

«» 0 THEN ERROR PCIB.BASERR 


ZERO.OHMS, SET.SPEED, J, J) 

<> 0 THEN ERROR PCIB.BASERR 

DCVOLTS, ACVOLTS, OHMS, R200MILLI) 

<> 0 THEN ERROR PCIB. BASERR 

R2, R20, R200, R2KILO) 

<> 0 THEN ERROR PCIB.BASERR 

R20KILO, R200KILO, R2MEGA, R2OMEGA) 

<> 0 THEN ERROR PCIB.BASERR 

AUTOM, R2.5, R12.5, J) 

<> 0 THEN ERROR PCIB. BASERR 
SET.COMPLEMENT, 


SBSE- DRIVER, OUTPUE NO.VALT, 


ENABLE.HANDSHAKE) 
GENE PETB.ERR <> O THEN ERROR PCIB.BASERR 
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CALL 


TI TART. BIT, 


68 
69 
“0 
71 
72 


84 


IF PCIB.ERR 
SME L.V(L, 
IF PCIB.ERR 
emer I.C(L, 
IF PCIB.ERR 
7 ı.c(I, 
IF PCIB.ERR 
EXER T.C(I, 
IF PCIB.ERR 
BAND I.C(I, 
IF PCIB.ERR 
$ 11 1.C(I, 
IF PCIB.ERR 
SME O I.C(L, 
IF PCIB.ERR 
I = 6 

G IT I.V(I 


FIT SYNMMETRY) 
55315: ECIB.ERR <> O THEN ERROR PCIB.BASERR 
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SALL I.V(I, 


I.V(I, 


, SET.FREQUENCY, 


DISABLE.HANDSHAKE, 
SET.NUM.BITS) 

«» 0 THEN ERROR PCIB.BASERR 
SET.LOGIC.SENSE, J, J, J) 
«» 0 THEN ERROR PCIB.BASERR 
POSITIVE, NEGATIVE, TWOS, UNSIGNED) 
«» 0 THEN ERROR PCIB.BASERR 

OC, TTL, RO, R1) 

«» 0 THEN ERROR PCIB.BASERR 

R2, R3, R4, R5) 

«» 0 THEN ERROR PCIB.BASERR 

R6, R7, R8, R9) 

«» 0 THEN ERROR PCIB.BASERR 

pelo sli, ri EIE, 

<> 0 THEN ERROR PCIB.BASERR 

R14, R15, R16, J) 

<> 0 THEN ERROR PCIB. BASERR 


SET. THRESHOLD, 


SEC AMPLITUDE; SET.OFFSET, 


SET.BURST.COUNT, 9, 9, 0) 


> 


87 
88 
89 
90 
2: 
92 
93 


IF PCIB.ERR 
CALL ul 
IF PCIB.ERR 
CALL tC em 
IF PCIB.ERR 


7 
CALL 


I. V (ür, 


<> O THEN ERROR PCIB.BASERR 

SINE, SQUARE, TRIANGLE, CONTINUOUS) 
<> O THEN ERROR PCIB.BASERR 

GATED, BURST, J, J) 

<> 0 THEN ERROR PCIB.BASERR 
AUTOSCALE, 


CALIBRATE, SET. SENSTIT ZI Z 


SETS VERI OFS) 


94 IF 


Per ELE 


95 CALL DV 
SET.DELAY) 
96 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 


SM 


CALL I 


<> 0 THEMERROR PCIB.BASERR 


SET.COUPLING, SET.POLARITY, SET.SWEEBSBEREE 


70 SET.TRIG.SOURCE, SET.TRIG.SLOPE, 


SET.TRIG.LEVEL, SET.TRIG.MODE) 
98 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 


o9 CAIL IE WIE, 


GET.SINGLE.WF, GET.TWO.WF, GET .VERT ih- 


GET.TIMEBASE. INFO) 
100 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 


101 


CAR, V ERM 


CALC.NF.STATS) 
102 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 


10:55 15 a” 


CALC. FREQUENCY ) 
104 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 


105 


CALL 


GET.TRIG. INFO, CALC.WFVOLT, CALE. WEIER 


CALC.RISETIME, CALC.FALLTIME, CALC£EE7777 


Iı. a CALC. PLUSWIDTH, CALC.MINUSWIDTH, 


CALG OVERSHOOT |) CALC PRESTOL) 
106 IF PCIB.ERR <> 0 THEN ERROR PEILBSDBASERR 


LO? 


SAL r, 


CALC.PK.TO.PK, SET.TIMEOUT, SCOPES a. 


MEASURE.SINGLE.WF) 


108 
OS 
110 
111 
PEZ 
ki 
114 
LES 
P Ro 
177 
116 
LL” 
120 
121 
122 
123 
124 
125 
126 
1277 


IF PCIB.ERR 
GAH ın 
IF PCIB.ERR 
CALL LA 
IF PCIB.ERR 
CALL I- CI 
IF PCIB.ERR 
CALL N GE 
IF PCIB.ERR 
CALL I C(I; 
IF PCIB.ERR 
GAME C(I; 
IF PCIB.ERR 
CALLA e 
IF PCIB.ERR 
CALL. ae 
IF PCIB.ERR 
CALL FEM 
IF PCIB.ERR 
ENAKE C(I, 


<> 0 THEN ERROR PCIB.BASERR 
MEASURE.TWO.WF, J, J, J) 

<> 0 THEN ERROR PCIB. BASERR 

RIONANO, R100NANO, RIMICRO, RIOMICRO) 
<> 0 THEN ERROR PCIB. BASERR 
R100MICRO, RIMILLI, R1OMILLI, 
<> 0 THEN ERROR PCIB. BASERR 
R1, R10, R20NANO, R200NANO) 
«» 0 THEN ERROR PCIB.BASERR 
R2MICRO, R20MICRO, R200MICRO, R2MILLI) 
«» 0 THEN ERROR PCIB.BASERR 

R20MILLI, R200MILLI, R2, R20) 

<> 0 THEN ERROR PCIB.BASERR 

R50NANO, R500NANO, R5MICRO, RSOMICRO) 
«» 0 THEN ERROR PCIB.BASERR 


R100MILLI) 


RSOOMICRO, RSMILLI, RSOMILLI, RSOOMILLI) 
<> © THEN ERROR PCIB. BASERR 

R5, R50, CHAN.A, CHAN.B) 

<> O THEN ERROR PCIB. BASERR 

EXTERNAL, POSITIVE, NEGATIVE, AC) 
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128 
129 
130 
131 
132 
133 
134 
135 
136 
137 
138 
199 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
800 
801 
899 
900 
901 
902 
903 
904 
iə 

205 
906 


IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

25015 17 C(L, DC, TRIGGERED, AUTO.TRIG, AUTO.LEVEL) 
IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

CALL I.C(I, X1, X10, STANDARD, AVERAGE) 

IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

I = 8 

CALL I.V(I, OPEN.CHANNEL, CLOSE.CHANNEL, J, J) 
IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

Ar rL C.S 

IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

IS = PCIB.DIR$ + "\PCIB.PLD" 

CALL L.P(IS) 

IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

IS = "DMM.01": I = 3: J = 0: K= 0: L= 1 

CALL DEFINE(DMM.O01, IS, I, J, K, L) 

IF PCIB.ERR <> O THEN ERROR PCIB. BASERR 

TS = "Func.Gen.01": I = 6: J = 0: K = 1: L= 1 
CALL DEFINE(Func.Gen.01, IŞ, I, J, K, L) 

IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

IS = "Scope.01": I = 7: J = O: K = 2: L= 1 
CALL DEFINE(Scope.01, I$, I, J, K, L) 

IF PCIB.ERR <> O THEN ERROR PCIB. BASERR 

£ 7 "Counter.01": I = 1: J = 0: K= 3: L= 1 

2 EE DEFINE(Counter.01, IŞ, I, J, K, L) 

IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

IS = "Dig.In.01": Is 4: Jz 0: K= 4: L= 1 
57775 DEFINE(Dig.In.01, I$, I, J, K, L) 

IF PCIB.ERR <> O THEN ERROR PCIB. BASERR 

mü “"Dig.Out.01": 1 4: J = 1: K = 4: L= 1 
CALL DEFINE(Dig.Out.01, I$, I, J, K, L) 

IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

IS = "Relay.Act.01": I = 8: J = 0: K = 5: L= 1 
WWW DEFINE(RELAY.ACT.01, I$, I, J, K, L) 

IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

IS = "Relay.Mux.01": I = 2: J = 0: K= 6: L= 1 
CALL DEFINE(RELAY.MUX.01, IS, I, J, K, L) 

IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

I$ = ENVIRONS("PANELS") + "\PANELS.EXE" 

ETE L.S(I$S) 

GOTO 2 

IF ERR = PCIB.BASERR THEN GOTO 903 

PRINT "BASIC error i"; ERR; " occurred in line "Ç; ERL 
STOP 

TMPERR = PCIB.ERR: IF TMPERR = O THEN TMPERR = PCIB.GLBERR 
PRINT "PC Instrument error #"; TMPERR; " detected at line 


ERL 


PRINT "Error: "; PCIB.ERR$ 
IF  LEFT$(PCIB.NAME$, 1)  <>  CHR$(32) THEN PRINT 


"Instrument: "; PCIB.NAMES 
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907 STOP 

908 COMMON PCIB.DIR$, PCIB.SEG 

909 COMMON LD.FILE, GET.MEM, PANELS, DEF.ERR 

910 COMMON PCIB.BASERR, PCIB.ERR, PCIB.ERR$, PCIB.NAMEŞ, 
PCIB.GLBERR 

9 ih il € O M M O N 
READ. REGISTER, READ. SELFID, DEFINE, INITIALIZE.SYSTEM, ENABLE. SY 
STEM, DISABLE.SYSTEM, INITIALIZE, POWER.ON,MEASURE, OUTPUT, START 
,HALT, ENABLE. INT. TRIGGER, DISABLE. INT. TRIGGER, ENABLE. OUTPUT, D 
ISABLE. OUTPUT, CHECK. DONE, GET. STATUS 

912 COMMON SET.FUNCTION, SET.RANGE, SET.MODE, WRITE.CAL, 
READ.CAL, STORE.CAL, DELAY,  SAVE.SYSTEM,  SET.GATETIME, 
SET.SAMPLES, SET.SLOPE, SET.SOURCE, ZERO.OHMS, SET.SPEED, 
SET.COMPLEMENT, SET.DRIVER, OUTPUT.NO.WAIT, ENABLE.HANDSHAKE, 
DISABLE. HANDSHAKE 

913 COMMON SET.THRESHOLD, SET.START.BIT, SET.NUM.BITS, 
SET.LOGIC.SENSE, SET.FREQUENCY, SET.AMPLITUDE, SET.OFFSET, 
SET.SYMMETRY, SET. BURST.COUNT, AUTOSCALE, CALIBRATE, 
SET.SENSITIVITY, SET.VERT.OFFSET, SET.COUPLING, SET.POLARITY, 
SET. SWEEPS PEED 

914 COMMON SET.DELAY, SET.TRIG.SOURCE, SET.TRIG.SLOPE, 
SET.TRIG.LEVEL, SET.TRIG.MODE, GET.SINGLE.WF, GET.TWO.WF, 
GET.VERT.INFO, GET.TIMEBASE.INFO, GET.TRIG.INFO, CALC.WFVOLT, 
CALC.WFTIME, CALC.WF.STATS, CALC.RISETIME, CALC.FALLTIME, 
CALC. PERIOD 

915 COMMON CALC. FREQUENCY, CALC.PLUSWIDTH, CALC.MINUSWIDTH, 
CALC.OVERSHOOT, CALC.PRESHOOT, CALC.PK.TO.PK, SET.TIMEOUT, 
SCOPE.START, MEASURE.SINGLE.WF, MEASURE.TWO.WF, OPEN.CHANNEL, 
CLOSE. CHANNEL 

916 COMMON FREQUENCY, AUTO.FREQ, PERIOD, AUTO.PER, INTERVAL, 
RATIO, TOTALIZE, R100MILLI, R1, R10, R100, R1KILO, R10MEGA, 
R100MEGA, CHAN.A, CHAN.B, POSITIVE, NEGATIVE, COMN, SEPARATE, 
DCVOLTS, ACVOLTS, OHMS, R200MILLI, R2, R20, R200, R2KILO, 
R20KILO _ 

, R200KILO 

917 COMMON R2MEGA, R20MEGA, AUTOM, R2.5, R12.5, POSITIVE, 
NEGATIVE, TWOS, UNSIGNED, OC, TTL, RO, R1, R2, R3, R4, R5, R6, 
R7, R8, R9, R10, R11, R12, R13, R14, R15, R16, SINE, SQUARE, 
TRIANGLE, CONTINUOUS, GATED,  BURST,  R10NANO,  R100NANO, 
RIMICRO, — 

R10MICRO, R100MICRO 

918 COMMON R1MILLI, R10MILLI, R100MILLI, R1, R10, R20NANO, 
R200NANO, R2MICRO, R20MICRO, R200MICRO, R2MILLI, R20MILLI, 
R200MILLI, R2, R20, RSONANO, R500NANO, RSMICRO, RSOMICRO, 
R500MICRO, RSMILLI, R50MILLI, R500MILLI, R5, R50, CHAN.A, 
CHAN.B, — 

EXTERNAL, POSITIVE 

919 COMMON NEGATIVE, AC, DC, TRIGGERED, AUTO.TRIG, AUTO.LEVEL, 
X1, X10, STANDARD, AVERAGE 
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920 COMMON DMM. 01, Func.Gen:01, 
Bin 01, Dig.Out.01, RELAV.ACT.O1, 
7 7 End PCIB Program Shell 

1000 REM This step initialzes the HP system 
FOTO CLS 


mora OPTION BASE 1 


Scope.01, 
RELAY.MUX.01 


Counter.01, 


15 7711 P(5), PA(50, 5), PP(50, 5), XPT(50), YPT(50), X(50), 
Y(50), YAW(50) 

1040 REM 

1050 CALL INITIALIZE.SYSTEM(PGMSHEL. HPC) 

1060 REM 

1070 REM SET FUNCTIONON THE 'DMM/! , "RELAY MUX , 'RELAY 
ACTUATOR” 

1080 REM 

1090 CALL SET.FUNCTION(DMM.01, DCVOLTS) 

1100 CALL SET.RANGE(DMM.01, AUTOM) 

1110 CALL DISABLE.INT.TRIGGER(DMM.01) 

1120 CALL ENABLE.OUTPUT(RELAY.MUX.01) 

1130 CALL ENABLE.OUTPUT(RELAY.ACT.01) 

1140 REM *X*XXXXXXXXXXX%X PROGRAM TRAVERSE XXKKAKKAXKKKAKKAKKAKKAKKAKAKAKAAKA 
1150 REM 

1160 REM OPEN THE COM PORT AND INITIALIZE THE MOTOR 
SETTINGS 

1170 OPEN "coml:1200,n,8,1,rs,cs,ds,cd" FOR RANDOM AS #1 
1180 REM SET MOTOR DEFAULT VALUES 

MED TA 2000, 2000,2000,2,2,2,0.000125,0.000125,0.000125 
£ UME RED V1, V2, V3, R1, R2, R3, C1, C2, C3 

1210 REM DEFINE CHARACTERS FOR DATA REDUCTION ALGORITHM 

1220 RN2Ş — "RENAME A:RAV.DAT " 

1230 HEADIŞ = " # X Y PI P2 P3 

P4 P5 YAW " 

1240 FORMATS = "787 ##-## ##-##  ###-###  ###.### FT. FHF 
HEH. HE #44 #44 #44. #4" 

1250 PRINT 

1260 PRINT UMkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
1270 PRINT "** USER MUST SELECT ’CAPS LOCK’ FUNCTION **" 
1280 PRINT tkkkxkkxkkkkkxkkkkkxkxkkkkkkkxkkkkkkxkkkkxkkkkkkkkkn 

1 “0 REM DISPLAY MOTOR DEFAULT SETTINGS 

1300 PRINT " xXkkkkikxkkkkkkkkkkkkkkkkkkkkkkkkkkokokoxk !l 
1310 PRINT " INITIALIZED VALUES FOR ALL MOTOR 
SETTINGS:' 

MIŻOSPRINT " VELOCITY — 1000 STEPS/SEC” 

1330 PRINT " RAMP (MOTOR ACCELERATION) = 2 (6000 
STEPS/SEC”2)" 

1340 PRINT " DEFAULT INCREMENTAL UNITS ARE 
INCHES" 

1350 PRINT " kkkkkkkkkikkkkkkkkkkkkkkkkkkkkkokkokkokin 
1360 PRINT 

1370 PRINT "NOTE!! USE MANUAL CONTROL TO INITIALIZE PROBE 
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POSTTION BEFORE” 


1380 
" 


1390 


PRINT >! SELECTING COMPUTER CONTROLLED MOVEMENT. 


PRINT 


1400 INPUT "MANUAL CONTROL OR COMPUTER CONTROL (ENTER “MAN” or 


“OP: 


1410 
1420 


CONS 
IF CONS = "CP" THEN 3490 
REM OPTION TO CHANGE DEFAULT SETTINGS OF VELOCITY OR 


ACCELERATION RAMP 


1430 
1440 
1450 


PRINT 
PRINT 
PRINTE: 


DO YOU WANT TO CHANGE THE VELOCITY OR 


ACCELERATION RAMP" 


1460 
1470 
1480 
THES 
1490 
1500 


PRINT " 
PRINT 
PRINT 


DEFAULT SETTINGS? (Y or N)" 


"IF ’NO’, THIS PROGRAM WILL THEN LET YOU DEFINE 
PRINT "DISTANCE YOU WANT TO MOVE (IN INCHES?. 157785 77 
PRINT "YOU CAN CHANGE ANY OR ALL OF THE DEFAULT SETTINGS 


FOR ANY MOTOR." 


1510 
1520 
1550 


PRINT 
PRINT 
PRINT 


1540 INPUT "DO YOU WANT TO CHANGE ANY OF THE DEFAULT SETTINGS? 


(Y "OE ST 


1550 
1560 
1570 
1580 
1590 
1600 
1610 


DS 
IF DS — ny 
IF DS = UN" 
REM 
REM 
PRINT 
PRINT 
INPUT 


THEN 1590 
UHEN 2220 


**** OPERATOR SELECTED MOTOR VARIABLES ****%* 


"WHICH DEFAULT VALUE? (ENTER '1'FOR VELOC OR “2” 


FOR ACCEL RAMP)"; L 


1620 
1630 
N)" 

1640 
1650 
1660 


ON L GOTO 1690, 1930 
PRINT "DO YOU WANT TO CHANGE THE DEFAULT VELOCITY? 


(Y OR 


INPUT V$ 
IF V$ = "Y" THEN 1690 
PRINT "DO YOU WANT TO CHANGE THE DEFAULT ACCELERATION 


RAMP? (Y or N)" 


1670 
1680 
1690 
1700 
1710 


"yu 
"NU 


THEN 1990 


IF R$ = 
= THEN 1450 


IF RS 
PRINT 
PRINT 
INPUT "WHICH MOTOR VELOCITY DO YOU WISH TO CHANGE? (1,2, 


OT s eg 


1720 
1730 
1740 


ON J GOTO 1730, 
PRINT 
PRINT 


18307212209 
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1750 
1760 
1770 
1780 
1790 


INPUT "ENTER DESIRED VELOCITY OF MOTOR #1"; V1 
PRINT 
PRINT 
PRINT 


PRINT "DO YOU WANT TO CHANGE VELOCITY OF ANOTHER MOTOR? 


(Y OR N)" 


1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
2910 
1920 
1930 
1940 
1950 
2, 
1960 
1970 
1980 
1990 
2000 
2010 
2020 


INPUT VŞ 
IF VS = 
IF VS = 
PRINT 
PRINT 
INPUT 
PRINT 
GOTO 1780 
PRINT 
PRINT 
PN PUT 
PRINT 
GOTO 1780 

PRINT 

PRINT 

INPUT "WHICH MOTOR ACCEL RAMP DO YOU WANT TO CHANGE? (1, 


"yn 
UN" 


THEN 1690 
THEN 1430 


TENTER DESIRED VELOCITY OF MOTOR 2"; V2 


TENTER T DESIRED VELOCITY OF MOTOR #3"; V3 


or 3)"; K 


ON K GOTO 1970, 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
PRINT 


2060, 2120 


"ENTER DESIRED ACCELERATION RAMP OF MOTOR #1"; R1 


"DO YOU WANT TO CHANGE THE ACCEL RAMP OF ANOTHER 


MOTOR (Y or N)?“ 


2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 


INPUT RMS 

IF RMŞ = "yt 
IF RMS = "N" 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
GOTO 2000 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
GOTO 2000 

REM 

REM DEFINE DISTANCE TO MOVE MOTOR 
PRINT 


THEN 1930 
THEN 1450 


"ENTER DESIRED ACCELERATION RAMP OF MOTOR #2"; R2 


"ENTER DESIRED ACCELERATION RAMP OF MOTOR #3"; R3 
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2210 PRINT 
2220 PRINT 
2230 REM INITIALIZE MOTOR INCREMENTS TO ZERO 
2240 I1 = 0 
2250 12 = 0 
2260 13 — 0 


2270 PRINT 
2 2 8 O Ë R I N T P 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
2290 PRINT " xx DEFINE WHICH MOTOR YOU WANT TO MOVE 

k k 
2300 PRINT " ** 

k ki 
2310 PRINT "xx NOTE!!! A POSITIVE (’+’) INCREMENT TO 
A MOTOR | xx" 
2320 PRINT " xx MOVES TRAVERSER AWAY FROM THAT PARTICULAR 
MOTOR kk 
2330 PRINT " xx 

de AL 
2340 PRINT " kk -- MOTOR #1 MOVES THE PROBE UPSTREAM AGAINST 
THE FLOW xx" 
2350 PRINT " ** -- MOTOR #2 MOVES THE PROBE TOWARD THE ACCESS 
WINDOW Ski 
2360 PRINT ha -- MOTOR #3 MOVES THE PROBE VERTICALLY 
DOWNWARD x 
2 3 7 0 P R I N T YT 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkii 
2380" PRENE 


2390 PRINT 

2400 INPUT "WHICH MOTOR DO YOU WANT TO MOVE? (1,2, OD MM 
2410 ON L GOTO 2420, 2680, 2970 

2420 PRINT 

2430 PRINT 

2440 PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #1?" 

2450 PRINT " kkkkkkkkk (ENTER DISTANCE IN INCHES) ********x*W" 
2460 INPUT Il 


2470 PRINT 

2480 PRINT W kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkka! 

2490 PRINT 

2500 PRINT "SUMMARY OF OPERATOR INPUTS:" 

2510 PRINT " MOTOR #1 VELOCITY = "; V1 

2520 PRINT " ACCELERATION RAMP = "; R1 
2530 PRENE INCREMENTAL DISTANCE = "> TT; 
INCHES? 

2540 PRINT " kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkki 

2550 PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y Or 
N)" 


2560 PRINT 
2570 PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER 1005 


66 


RETURN" 


2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2720 
2740 
2750 
2760 
2770 
2780 
2/90 
2800 


PRINT 
INPUT 
TE VS 
GOSUB 
PRINT 
PRINT 
INPUT 
TE CS 
EE CS 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 


"TO VARIABLE SELECTION SUBROUTINE." 
VS 

= "ун 
3410 


THEN 1430 


"DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N)?" 
ES 

= "ун 
z UN" 


THEN 2220 
THEN 3260 


"HOW FAR DO YOU WANT TO MOVE MOTOR #2?" 
" xkkkkkkk* (ENTER DISTANCE IN INCHES) *****X*XXX" 
bo 


REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


11: HES" 


2810 
2820 
2530 
2840 
N)" 

2850 
2860 


PRINT 
PRINT 
PRINT 
PRINT 


PRINT 
PRINT 


RETURN" 


2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 


PRINT 
INPUT 
m VS 
GOSUB 
PRINT 
PRINT 
INPUT 
IR CS 
HE (CS 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 


n kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk! 


"SUMMARY OF OPERATOR INPUTS:" 


Š MOTOR #2 VELOCITY - "; V2 
2 ACCELERATION RAMP = "; R2 
E INCREMENTAL DISTANCE - "; I2; 


O kxkxkxkxkxkxkxkkkkkkkkkkkkkkkkkkkkkkkkk 


"DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 


CENTER 11107 “TEL MOTOR MOVEMENT. CENTER Y” TO 


"TO VARIABLE SELECTION SUBROUTINE. " 
VS 


= "yu 
3410 


THEN 1430 


"DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? mor Mu! 
CS 

= "ун 
= "Ы" 


THEN 2220 
THEN 3260 


"HOW FAR DO YOU WANT TO MOVE MOTOR #3?" 
"ökkkkkkkkk (ENTER DISTANCE IN INCHES) 
13 


kkkkkkkkk" 
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3050 
3040 
3050 
3060 
3070 
3080 
3090 


REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


“INCHES: 


3100 
3110 
3120 
3130 
3140 
N)" 

3150 
3160 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


PRINT 
PRINT 


RETURN" 


3170 
3180 
3190 
3200 
3210 
3220 
3230 


PRINT 
INPUT 
IF VS 
GOSUB 
PRINT 
PRINT 
INPUT 


Cy won" 


3240 


IF MS 


ENTER 


3250 PRINT 
3260 PRINT 
MOVEMENT?" 
3270 PRINT 
3280 PRINT 
ZEROLZED." 
3290 PRINT 
MOVEMENT." 
3300 PRINT 
3310 PRINT 
3320 INPUT 
Or N) “ANS 


3330 
3340 
5350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 


":C3="; 


IF N$ 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
END 

REM 

PRINT 
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ü kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk! 


"SUMMARY OF OPERATOR INPUTS:" 


m MOTOR #3 VELOCITY = "> V3 
7 ACCELERATION RAMP = ", R3 
` INCREMENTAL DISTANCE = "; I3; 


N kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk" 


"DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or 


”N” TO START MOTOR MOVEMENT. ENTERS IS 


"TO VARIABLE SELECTION SUBROUTINE." 
VS 


nyn 
3410 


THEN 1430 


"DO YOU WANT TO INPUT ANOTHER MANUAL MOTOR MOVEMENT 
MS 


"VU RENER 210 


"DO YOU WANT TO INPUT COMPUTER CONTROLLED MOTOR 


" *kkkkkkkk NOTE!!! kkkkkkkkk mn 
" ALL PREVIOUS MOTOR INCREMENT INPUTS HAVE BEEN 


"PROGAM WILL LET YOU CHOOSE MANUAL OR CP-CONTROLLED 


"kkkkk (IF kök kk kü 


‘NO’, THE PROGRAM WILL END). 


"DO YOU WANT COMPUTER CONTROLLED MOTOR MOVEMENT (Y 


aya THEN 3500 


" kkkkkkkkkkkkkkkkkkkkkkkkkkkk! 


Ы THE PROGRAM HAS ENDED." 
" * * * he ke  * * * k k k k k k kk k kk kk k i 


******* MOTOR MOVEMENT SUBROUTINE ********* 


Yi, "8": PRINT #1, "E"; "c1"; C1; HG 'U R 
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1255 275 501 ?1, "Ile", Il; ":Vl1z"; V1; ":Ri="; Ri; 


EER EE TINT $1, ":I2-"; 12, ":V2-", V2; ":R2="; R2 

RE OMERINT #1, “T3="; I3; ":V3="; V3; ":R3="; R3; ":@" 
3460 RETURN 

3470 REM kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkxkk 
3480 REM kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkxk 
3490 PRINT 

3500 REM ******* COMPUTER CONTROLLED MOVEMENT *******x 
3510 PRINT 


3520 PRINT "THE PRESSURE DATA WILL BE WRITTEN TO FILES ON 
DRIVE “A? " 

23530 PRINT 

3540 PRINT "YOU WILL BE ASKED TO INPUT FILE NAMES FOR THESE." 
55:22) PRINT 

SES OM NPUT "IS A FORMATTED DISK IN DRIVE “A”? PRESS “ENTER” 
MONCONTINUE"; DS 


2570 PRINT 
23530 PRINT 
520 PRINT 
3600 PRINT " * de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de e de de Re NM 
BONOMPRINT ' kök NOBEN! kk! 
OZO" PRINT " kk COMPUTER CONTROLLED MOVEMENT kk! 
BOJSONPRINT " kok IS PROGRAMMED WITH A ern 
3640 PRINT " ** DEFAULTED NEGATIVE MOTOR INCREMENT xx" 
26550 PRINT " ** (i.e. MOTOR £3 WILL MOVE UPWARD kk! 
3660 PRINT " kk BY ENTERING A (+) DISTANCE). Ski 
3670 PRINT " de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de Ae de PI 
3680 PRINT 


3690 REM SET INITIAL MOVEMENT DISTANCE AND NUMBER OF DATA 
POINTS TO ZERO 

3700 HT = 0 

3710 WD = 0 

3720 DIST = O 

3730 XPT = 0 

3740 YPT = 0 

3750 N = 0 

3760 PRINT 

3770 PRINT 

3780 INPUT "WHAT IS THE DIMENSION ( X , Y ) (IN INCHES) THAT 
YOU WANT TO MEASURE."; WD, HT 

3790 PRINT 

3800 INPUT "WHAT IS THE STEP (IN INCHES) THAT YOU WANT TO 
MOVE Hy DIST 

1 ET = INT(HT / DIST) 
507 ET — INT(WD / DIST) 
3830 N — XPT * YPT 

3840 PRINT 

SPRINT “THERE ARE "; XPT; " * "; YPT; "Ç" = "=+ N; " POINTS 
TO BE MEASURED " 


dS xd 
+ 1 
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3860 
3870 
3880 
3890 
3900 
3910 
3920 


NAMES 


3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 


EXTENSION) "; 


4040 
4050 
4060 
4070 
4080 


PRINT | 

INPUT "ARE THE NUMBER OF POINTS IS OK.(Y OR N)"; CS 

IF C$ = "N" THEN 3780 

CLS 

N = XPT 

IF (N < 1) OR (N > 99) GOTO 3780 

REM *** GENERATING STRING STRING SEGMENTS FOR DATA FILE 


REM ** STRING NUMBER "1" 
REM ** ENDING STRING NUMBER 


BŞ MIDS(STRS(1), 
ES MIDŞ (STRŞ (N), 
XS OOO 

EXS =m. DAP! 

CLS 
PRINT 
PRINT 
PRINT (XS BS F EXS) "U 
PRINT 
PRINT 
INPUT 


2): 


2) s "N" 


"DATA FILES WILL BE INCREMENTED FROM:" 


To "MORS EES ad 


"ENTER NO 


F2$ 


DATA FILE NAME (6 CHARACTERS MAX -- 
PRINT 

PRINT 

IF LEN(F2$) > 6 OR LEN(F2$) < 1 GOTO 4030 

CLS 


PRINT N; "DATA FILES WILL BE GENERATED AND INCREMENTED AS 


FOLLOWS: " 


4090 
4100 
4110 
4120 
4130 
4140 
N) n 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT "ARE THE NUMBER OF POINTS AND FILE NAMES OK.(Y OR 
ÉS 

IF C$ = "N" GOTO 3780 

IF C$ - "Y" GOTO 4180 

GOTO 4140 

CLS 

PRINT 

PRINT 

REM SET INITIAL POSITION DATA 

X(1) = -DIST 

ma DIST 

FOR IX = 2 TO XPT +1 

X(IX) = 0 
NEXT IX 
FOR JY = 
Y(IY) = 0 
NEXT JY 
FOR I = 


(F2$ + BS + EXS); " To m": (F2S + ES (967 


Z EO Pi MM 


0) MET 
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52502 11 
4304 12 
4306 13 
n- -HOR J = 1 TO YPT 

4320 REM MOTOR CP-CONTROLLED MOTOR MOVEMENT 


O 
O 
O 


2230 T1 = О 

2570012 -0 

4350 I3 = O 

4360 REM EACH POINT TAKE 10 TIMES READINGS 
4370 X(I + 1) = X(I) + DIST 

4380 XPT(J) = X(I + 1) 

4390 Y(J + 1) = Y(J) + DIST 


4400 YPT(J) = Y(J + 1) 

4405 INPUT " ADJUST THE WHEEL TO MAKE THE P2 =P3,INPUT THE YAW 
ANGLE"; YAW(J) 

4408 PRINT 

4410 INPUT " PRESS “ENTER” TO START THE MEASUREMENT"; MOVES 

4420 REM 

4430 REM READ FIVE CHANNELS AND DISPLAY THE DATA 

4440 REM 

4450 STEPPER = 4 

4460 SWITCH = 3 

4470 HOMER = 8 

2190 DELAYI = .1 
4490 DELAY2 = 1 

4500 REM SET THE S.V PORT TO #4 

4510 FOR IL = 1 TO 3 

4520 THYME = TIMER 

4530 CALL OUTPUT (RELAY.ACT.01,STEPPER) 

4540 CHKTIME = TIMER 

4550 IF CHKTIME < (THYME + DELAY1) GOTO 4540 
4560 CALL OPEN.CHANNEL(RELAY.ACT.01, SWITCH) 
4570 CLS 

4580 NEXT IL 

4590 PRINT 

4600 PRINT " NOW IS POINT "; J 

4610 REM START MEASURE FROM PORT 4 TO PORT 8 
4620 FOR JJ = 1 TO 5 

4630 CALL OUTPUT(RELAY.ACT.01,STEPPER) 

4640 CHKTIME = TIMER 

4650 IF CHKTIME < (THYME + DELAY2) GOTO 4640 
4660 REM EACH PORT SAMPLE 10 TIMES 

4670 FOR II = 1 TO 10 

4680 ROUT = 1 

4690 CALL OUTPUT(RELAY.MUX.01,ROUT) 

4700 CALL MEASURE(DMM.01, VOLTS) 

4710 PA(II, JJ) = VOLTS 

4720 NEXT II 

4730 CALL OPEN.CHANNEL(RELAY.ACT.01, SWITCH) 


FAJL 


4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4 


IF JJ = 
NEXT JJ 
REM HOME THE S.V PORT TO #48 

CALL OUTPUT (RELAY.ACT.01,HOMER) 


5 THEN 4760 


CALL OPEN.CHANNEL(RELAY.ACT.01, HOMER) 

REM 

REM DISPLAY THE SAMPLE DATA 

REM 

PRINT HEADIŞ 

FOR IS= 1 TO 10 

8 4 O P R I N si U S I N G 


FORMATS ;IS,XPT(J),YPT(J) ,PA(IS,1),PA(IS,2),PA(1IS,3) ,PA(IS, 4) 
,PA(IS,5) ,YAW(J) 


4850 NEXT IS 

4860 REM 

4870 REM AVERAGE THE DATA 

4880 REM 

4890 FOR JA < 1 TO 5 

4900 TOTAL €70 

4910 FOR IA = 1 TO 10 

4920 TOTAL = TOTAL + PA(TA, JA) 

4930 NEXT IA 

4940 AVERAGE = TOTAL / 10 

4950 P(JA) = AVERAGE 

4960 NEXT JA 

4970 PRINT 

4980 PRINT "THE AVERAGES ARE: " 

5000 PRINT HEAD1$ 

5010 FOR JD = 1 TO 5 

5020 PP(J, JD) = P(JD) 

5030 NEXT JD 

5040 PRINT USING FORMATS; J; XPT(J); YPT(J); PP(J, 1); PP(J, 
2); PP(J, 3); PP(J, 4); PP(J, 5); YAW(J) 

5045 PRINT 

5048 PRINT USING "THE NULLING ERROR IS +#.####"; PP(J, 3) - 
PP(J, 2) 

5049 PRINT 

5050 PRINT "DO YOU WANT RE-MEASURE AGAIN (Y / N)" 
5060 PRINT 

5062 PRINT "IF “Y” WILL RE-SAMPLE AGAIN." 

5064 PRINT 

5070 INPUT "IF 'N' WILL MOVE THE TRAVERSER STEP UPWARD (WAIT 
7 SEE "CS 

5075 PRINT 

5080 IF C$ € "Y" THEN 4405 

5082 IF CS = "N" THEN 5090 

5084 GOTO 5070 

5090 IF J = YPT THEN 5160 

5100 REM 
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5110 


REM MOVE THE TRAVERSER STEP UPWARD. 


5120 REM 

5130 13 = -DIST 

5140 GOSUB 3410 

5150 NEXT J 

5160 REM*** STORE DATA BEFORE NEXT SAMPLE*** 

5170 OPEN "A:\RAW.DAT" FOR OUTPUT AS #2 

5180 PRINT #2, HEAD1$ 

5190 FOR ID = 1 TO YPT 

5200 PRINT #2, USING FORMATS; ID; XPT(ID); YPT(ID); PP(ID, 1); 
PP(ID, 2); PP(ID, 3); PP(ID, 4); PP(ID, 5); YAW(ID) 
5210 NEXT ID 

5220 CLOSE #2 

5230 REM *** GENERATING INCREMENTED DATA FILE NAME 

5240 IF (I > 10) OR (I = 10) THEN IS = MIDS(STRS(I), 2) 
5250 IF (I < 10) THEN IS = (MIDS(STRS(O), 2) + MIDS(STRS(I), 
2)) 

5260 FI2$ = (F2$ + IS + EX$) 

5270 PRINT 

5280 PRINT " WRITING DATA FILE "; FI2$ 

5290 DF2$ = RN2$ + FI2$ 

5300 REM ** RENAME DATA FILE 

5310 SHELL DF2$ 

5320 REM 

5330 REM MOVE THE TRAVERSER TO THE NEXT SAMPLE POSITION 


5340 REM 

5350 PRINT 

5360 IF I = XPT THEN 5430 

5370 INPUT "THEN PRESS ’ENTER’ FOR NEXT COLUMN SAMPLE( 90 SEC) 
MOVES 

39012 = -DIST 

5400 I3 = HT 

5410 GOSUB 3410 

5420 NEXT I 

5430 CLS 

5440 PRINT "ALL MOVEMENTS COMPLETE" 

5450 PRINT 

5460 PRINT 

5470 PRINT "YOU VANT TO REPOSITION TRAVERSER FOR ANOTHER 
MOVEMENT (Y OR N)?" 

5480 PRINT 

5490 PRINT "IF “Y”, THE PROGRAM VILL TAKE YOU TO MANUAL 
CONTROL SUBROUTINE." 

5500 PRINT "IF 'N', THE PROGRAM WILL END." 

5510 PRINT 

5520 INPUT "ANOTHER MOVEMENT"; RS 

5530 IF RS = "Y" THEN 1370 

577702 1F RŞ — "N" THEN 3370 
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APPENDIX B 


1 DEF SEG : CLEAR , &HFEOO: GOTO 4’Begin PCIB Program Shell 
2 GOTO 1000 ” User program 

3 GOTO 900 * Error handling 

4 I = &HFEOO’ Copyright Hewlett-Packard 1984,1985 
5 PCIB.DIRŞ = ENVIRONS ("PCIB") 

6 IS = PCIB.DIRS + "\PCIBILC.BLD" 

7 BLOAD IS, I 

8 CALL I(PCIB. DIRS, 12,272 PCB Hel er 

9 IF J% = O THEN GOTO 13 

10 PRINT “Unable to load."; 

11 PRINT " (Error 7”? 4:55 

12 END 

1 


14 DEF SEG = PCIB.SEG: O.S = 5: C.S = 10: I.V = 15 

15 I.C = 20: L.P = 25: LD.FILE = 30 

16 GET.MEM = 35: L.S = 40: PANELS = 45: DEF.ERR = 50 

17 PCIB.ERR$ = STRING$(64, 32): PCIB.NAME$ < STRING$(16, 32) 
18 CALL DEF.ERR(PCIB.ERR, PCIB.ERR$, PCIB.NAME$, PCIB.GLBERR): 
PCIB.BASERR = 255 

19 ON ERROR GOTO 3 

20 J = -1 

21 15 = PCIB.DIR$ + "\PCIB.SYN" 

INCA OS (ES) 

23 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

24 I = 0 

25 u CALL ии READ.REGISTER, READ.SELFID, DEFINE, 
INITIALIZE.SYSTEM) 

26 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

27 CALL I.V(I, ENABLE.SYSTEM, DISABLE.SYSTEM, INITIALIZE, 
POVER.ON) 

28 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

29 CALL I.V(I,MEASURE,OUTPUT,START,HALT) 

30 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

31 CALL I.V(I,  ENABLE.INT.TRIGGER,  DISABLE.INT.TRIGGER, 
ENABLE.OUTPUT, DISABLE.OUTPUT) 

32 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

33 CALL I.V(I,  CHECK.DONE, GET.STATUS,  SET.FUNCTION, 
SET.RANGE) 

34 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

35 CALL I.V(I, SET.MODE, WRITE.CAL, READ.CAL, STORE CA 

36 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

37 CALL I.V(I, DELAY, SAVE.SYSTEM, J, J) 

38 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

39 I = 1 

40 CALL I.V(I, SET.GATETIME, SET.SAMPLES, SET.SLOPE, 
SET.SOURCE) 
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63 
64 
65 


IF PCIB.ERR 
SMEL I.C(L, 
IF PCIB.ERR 
$ EE T.CÇI, 
IF PCIB.ERR 
S ER T.C(I, 
IF PCIB.ERR 
CALL I.C(I, 
IF PCIB.ERR 
SMEM I.C(I, 
IF PCIB.ERR 
I = 2 

I = 3 

© I.V(I, 
IF PCIB.ERR 
eon I.C(I, 
IF PCIB.ERR 
270 I.C(I, 
IF PCIB.ERR 
eM L.C(I, 
IF PCIB.ERR 
em I.C(I, 
IF PCIB.ERR 
I = 4 

£ EE I.V(I, 


<> 0 THEN ERROR PCIB.BASERR 
FREQUENCY, AUTO.FREQ, PERIOD, AUTO. PER) 
<> 0 THEN ERROR PCIB.BASERR 

INTERVAL, RATIO, TOTALIZE, R100MILLI) 
«» 0 THEN ERROR PCIB.BASERR 

R1, R10, R100, RIKILO) 

<> 0 THEN ERROR PCIB.BASERR 

R10MEGA, R100MEGA, CHAN.A, CHAN.B) 

<> 0 THEN ERROR PCIB.BASERR 

POSITIVE, NEGATIVE, COMN, SEPARATE) 

<> 0 THEN ERROR PCIB.BASERR 


ZERO.OHMS, SET.SPEED, U, J) 
<> O THEN ERROR PCIB.BASERR 


DCVOLTS, ACVOLTS, OHMS, R200MILLI) 
<> 0 THEN ERROR PCIB. BASERR 
GEET: 0200) FOETLO) 


<> 0 THEN ERROR PCIB. BASERR 
R20KILO, R200KILO, R2MEGA, R20MEGA) 
<> 0 THEN ERROR PCIB.BASERR 
AUTOM, R2.5, R12.5, J) 
<> 0 THEN ERROR PCIB. BASERR 

SET. COMPLEMENT, 


SET DEEVER, OUTPUT.NO.WAIT, 


ENABLE.HANDSHAKE) 
CONTE PCIB.ERR <> O THEN ERROR PCIB. BASERR 
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CALL 


SET.START.BIT, 


82 
83 
84 


IF PCIB.ERR 
F TIL TI.V(I, 
IF PCIB.ERR 
SN T.C(I, 
IF PCIB.ERR 
EXER T.C(I, 
IF PCIB.ERR 
2556 I.c(I, 


LE PEITB.ERR 


Beer, 
IF PCIB.ERR 
2 0 T. C(I, 
IF PCIB.ERR 
SIN lé 
IF PCIB.ERR 
I = 6 


“SEE I.V(LI, 


SET.SYMMETRY) 
sie alt ECIB.ERR <> 0 THEN ERROR PCIB.BASERR 
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BAEL I.V(I, 


Fr 


DISABLE.HANDSHAKE, 
SET.NUM.BITS) 

«» 0 THEN ERROR PCIB.BASERR 
SET.LOGIC.SENSE, J, J, J) 
«» 0 THEN ERROR PCIB.BASERR 
POSITIVE, NEGATIVE, TWOS, UNSIGNED) 
«» 0 THEN ERROR PCIB.BASERR 

OC, TTL, RO, R1) 

«» 0 THEN ERROR PCIB.BASERR 

R2, R3, R4, R5) 

«» 0 THEN ERROR PCIB.BASERR 

R6, R7, R8, R9) 

«» 0 THEN ERROR PCIB.BASERR 

R10, R11, R12, R13) 

«» 0 THEN ERROR PCIB.BASERR 

R14, R15, R16, J) 

<> O THEN ERROR PCIB. BASERR 


SET.THRESHOLD, 


SET PPEOUEPICY OR O AMPDITUDE, SEI-OFFSET, 


SET.BURST.COUNT, J, J, J) 


75 


93 


IF PCIB.ERR 
CALL ree (ae 
IF PCIB.ERR 
CALLS Ae 
IF PCIB.ERR 
I = 7 


CALL 


TV 


<> 0 THEN ERROR PCIB.BASERR 

SINE, SQUARE, TRIANGLE, CONTINUOUS) 
<> 0 THEN ERROR PCIB. BASERR 

GATED, BURST, J, J) 

<> 0 THEN ERROR PCIB. BASERR 
AUTOSCALE, 


CALIBRATE, SET .SENS TRAAN 


SET.VERT.OFFSET) 

94 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 
QONAR V (1! 
SET.DELAY) 

96 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 


dı 


CALL I 


SET.COUPLING, SET.POLARITY, SET.SWEEPSPEEDS 


757: SET.TRIG.SOURCE, SET.TRIG.SLOPE, 


SET. TRIG. LEVEL, SET SEIREIGAMODE, 
98 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 


39 


CALL IVO 


GET.SINGLE.WF, GET.TWO.WF, GET .VERTM AFR 


GET. TIMEBASE. THEO) 
100 IF PCIB.ERR €” 0 THEN ERROR PCIB.BASERR 


101 


CAGE, TE. VOE 


CALC.WF.STATS) 
102 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 


OS Er mn Vr. 


CALC.FREQUENCY) 
104 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 


105 


CALL 


GET.TRIG.INFO, CALC.WFVOLT, CALE WE RUM 


CALC.RISETIME, CALC.FALLTIME, CALC £177:72 


IV i; CALC.PLUSWIDTH, CALC.MINUSWIDTH, 


CALC.OVERSHOOT, CALC. PRESHOOT) 
106 IF PCIB.ERR <> OMEJEN ERROR, POMENE BER 


io EER IE Or 


CALC.PK.TO.PK, SET .TIMEOUT, SCOPE TARIE 


MEASURE . SINGLE .WF) 


108 
109 
110 
111 
172 
115 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 


IF PCIB.ERR 
CALM (E 
IF PCIB.ERR 
FALL: ISCHTE 
IF PCIB.ERR 
Fre 
IF PCIB.ERR 
елини GC C. 
IF PCIB.ERR 
CARIBE C (NDS 
IF PCIB.ERR 
CAM CI: 
IF PCIB.ERR 
CALL ME (di 
IF PCIB.ERR 
CALL T C 
IF PCIB.ERR 
CAMPER GM 
IF PCIB.ERR 
CALL EET 


<> 0 THEN ERROR PCIB.BASERR 
MEASURE.TWO.WF, J, J, J) 

<> O THEN ERROR PCIB.BASERR 

R10NANO, R100NANO, RIMICRO, R10MICRO) 
«» 0 THEN ERROR PCIB.BASERR 
RIOOMICRO, RIMILLI, RIOMILLI, 
<> 0 THEN ERROR PCIB.BASERR 
R1, R10, R20NANO, R200NANO) 
«» 0 THEN ERROR PCIB.BASERR 
R2MICRO, R20MICRO, R200MICRO, R2MILLI) 
«» 0 THEN ERROR PCIB.BASERR 

R20MILLI, R200MILLI, R2, R20) 

«» 0 THEN ERROR PCIB.BASERR 

R5ONANO, R500NANO, RSMICRO, R50MICRO) 
<> 0 THEN ERROR PCIB.BASERR 
R500MICRO, RSMILLI, RSOMILLI, 
<> 0 THEN ERROR PCIB.BASERR 
R5, R50, CHAN.A, CHAN.B) 

<> 0 THEN ERROR PCIB.BASERR 
EXTERNAL, POSITIVE, NEGATIVE, AC) 


R100MILLI) 


R500MILLI) 
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128 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

5 7:36 I.C(I, DC, TRIGGERED, AUTO.TRIG, AUTO. LEVEL) 
130 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

131 CALL I.C(I, X1, X10, STANDARD, AVERAGE) 

132 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

133 I = 8 

134 CALL I.V(I, OPEN.CHANNEL, CLOSE.CHANNEL, J, J) 
135 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

136 CALL C.S 

137 IF PCIB.ERR <> O THEN ERROR PCIB. BASERR 

138 IS = PCIB.DIRS + " E€CLB. PED" 

139 CALL L.P(IS) 

140 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 

141 IS = "DMM.01": I = 3: J = 0: K = 0: L = 1 

142 CALL DEFINE(DMM.01, IS, I, J, K, L) 

143 IF PCIB.ERR <> O THEN ERROR PCIB. BASERR 

144 IS = "Func.Gen.01": TI = 6: J = 0: K = 1: L = 1 
145 CALL DEFINE(FUNC.GEN.01, I$, I, J, K, L) 

146 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

147 I$ < "Scope.01": I < 7: J < 0: K < 2: L < 1 
148 CALL DEFINE(SCOPE.O1, I$, I, J, K, L) 

149 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

150 IS — "Counter.01": I = 1: J = 0: K= 3: L = 1 
151 CALL DEFINE(COUNTER.O1, I$, I, J, K, L) 

152 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 
ISS "Dig.In.01": I = 4: J= 0: K = 4: L= 1 
154 CALL DEFINE(DIG.IN.O1, I$, I, J, K, L) 

155 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

156 IŞ € "Dig.Out.O1": IT = 4: J = 1: K= 4: L= 1 
25.7 DEFINE(DIG.OUT.01, I$, I, J, K, L) 

158 IF PCIB.ERR <> 0 THEN ERROR PCIB. BASERR 


159 IS = "Relay.Act.01": I = 8: J = 0: K= 5: L= 1 

1777 :7£ DEFINE(RELAY.ACT.O1, IŞ, I, J, K, L) 

161 IF PCIB.ERR <> O THEN ERROR PCIB.BASERR 

162 IS < "Relay.Mux.01": I < 2: J = 0: K= 6: L= 1 

NSENMCHER DEFINE(RELAY.MUX.O1, IS, I, J, K, L) 

164 IF PCIB.ERR <> 0 THEN ERROR PCIB.BASERR 

800 IS = ENVIRONS("PANELS") + "\PANELS.EXE" 

220055 L.S(IS) 

899 GOTO 2 

900 IF ERR = PCIB.BASERR THEN GOTO 903 

5012059)? "BASIC error f"? ERR, " occurred in line "; ERL 
902 STOP 

903 TMPERR = PCIB.ERR: IF TMPERR = O THEN TMPERR = PCIB.GLBERR 
904 PRINT "PC Instrument error #"; TMPERR; " detected at line 
Ue ERL 

905 PRINT "Error: "; PCIB.ERRS 


906 IF  LEFT$(PCIB.NAME$, 1) <> CHR$(32) THEN PRINT 
"Instrument: "; PCIB. NAMES 


NU 


907 STOP 
908 COMMON PCIB.DIR$, PCIB.SEG 

909 COMMON LD.FILE, GET.MEM, PANELS, DEF.ERR 

910 COMMON PCIB.BASERR, PCIB.ERR, PCIB.ERRS, PCIB.NAMES, 
PCIB.GLBERR 

9 1 1 C O M M O N 
READ. REGISTER,READ.SELFID,DEFINE,INITIALIZE.SYSTEM,ENABLE.SY 
STEM, DISABLE.SYSTEM, INITIALIZE, POWER.ON, MEASURE, OUTPUT, START 
, HALT, ENABLE. INT. TRIGGER, DISABLE. INT. TRIGGER, ENABLE. OUTPUT, D 
ISABLE. OUTPUT, CHECK. DONE, GET. STATUS 

912 COMMON SET.FUNCTION, SET.RANGE, SET.MODE, WRITE.CAL, 
READ.CAL, STORE.CAL, DELAY, SAVE.SYSTEM, SET.GATETIME, 
SET.SAMPLES, SET.SLOPE, SET.SOURCE, ZERO.OHMS, SET.SPEED, 
SET.COMPLEMENT, SET.DRIVER, OUTPUT.NO.WAIT, ENABLE.HANDSHAKE, 
DISABLE .HANDSHAKE 

913 COMMON SET.THRESHOLD, SET.START.BIT, SET.NUM.BITS, 
SET.LOGIC.SENSE, SET.FREQUENCY, SET.AMPLITUDE, SET.OFFSET, 
SET.SYMMETRY, SET.BURST.COUNT, AUTOSCALE, CALIBRATE, 
SET.SENSITIVITY, SET.VERT.OFFSET, SET.COUPLING, SET.POLARITY, 
SET.SWEEPSPEED 

914 COMMON SET.DELAY, SET.TRIG.SOURCE, SET.TRIG.SLOPE, 
SET.TRIG. LEVEL, SET.TRIG.MODE, GET.SINGLE.WF,  GET.TWO.WF, 
GET.VERT.INFO, GET.TIMEBASE.INFO, GET.TRIG.INFO, CALC.WFVOLT, 
CALC.WFTIME, CALC.WF.STATS, CALC.RISETIME, CALC.FALLTIME, 
CALC. PERIOD 

915 COMMON CALC.FREQUENCY, CALC.PLUSWIDTH, CALC.MINUSWIDTH, 
CALC.OVERSHOOT, CALC.PRESHOOT, CALC.PK.TO.PK, SET.TIMEOUT, 
SCOPE.START, MEASURE.SINGLE.WF, MEASURE.TWO.WF, OPEN.CHANNEL, 
CLOSE. CHANNEL 

916 COMMON FREQUENCY, AUTO.FREQ, PERIOD, AUTO.PER, INTERVAL, 
RATIO, TOTALIZE, R100MILLI, R1, R10, R100, RIKILO, R10MEGA, 
R100MEGA, CHAN.A, CHAN.B, POSITIVE, NEGATIVE, COMN, SEPARATE, 
DCVOLTS, ACVOLTS, OHMS, R200MILLI, R2, R20, R200, R2KILO, 
R20KILO _ 

, R200KILO 

917 COMMON R2MEGA, R20MEGA, AUTOM, R2.5, R12.5, POSITIVE, 
NEGATIVE, TWOS, UNSIGNED, OC, TTL, RO, Rl, R2, R3, R4, R5, R6, 
R7, R8, R9, R10, R11, R12, R13, R14, R15, R16, SINE, SQUARE, 
TRIANGLE, CONTINUOUS, GATED,  BURST,  R10NANO,  R100NANO, 
RIMICRO, — 

R10MICRO, R100MICRO 

918 COMMON RIMILLI, RIOMILLI, RIOOMILLI, R1, R10, R20NANO, 
R200NANO, R2MICRO, R20MICRO, R200MICRO, R2MILLI, R20MILLI, 
R200MILLI, R2, R20, RSONANO, R500NANO, R5MICRO, R50MICRO, 
R500MICRO, R5MILLI, RSOMILLI, R500MILLI, RS, R50, CHAN.A, 
CHAN.B, _ 

EXTERNAL, POSITIVE 

919 COMMON NEGATIVE, AC, DC, TRIGGERED, AUTO.TRIG, AUTO.LEVEL, 
X1, X10, STANDARD, AVERAGE 
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920 COMMON DMM.01, FUNC GENS O1; SCOPE OL, COUNTER. 01, 
ик ит 01, DIG.OUT.O1, RELAY.ACT.O1, RELAY .MUX.O01 

99 End PCIB Program Shell 

1000 KEM 

1010 REM This step initialzes the HP system 

1020 CLS 

1030 OPTION BASE 1 

1040 DIM P(10), PA(50, 6), PP(50, 6), XPT(40), CAL(40) 
1050 CALL INITIALIZE.SYSTEM(PGMSHEL.HPC) 

1060 REM 

1070 REM All PC devices nov have an initial state 

1080 REM Set function on the DMM and Relay MUX 

1090 REM 

1100 CALL SET.FUNCTION(DMM.01, DCVOLTS) 

1110 CALL SET.RANGE(DMM.01, AUTOM) 

1120 CALL DISABLE. INT.TRIGGER(DMM.01) 

1130 CALL ENABLE.OUTPUT(RELAY.MUX.01) 

1140 FORMATS = "## tt. HH Hit RR. RARR gn, RTH Ht. Hitt i 
HE- PREE FF. FTE" 

1200 FOR I = 1 TO 10 

1210 CAL(I) = 0! 

220 NEXT I 

1510 REM 

1520 REM READ THE VOLTAGE OF 48TH CHANNEL AND DISPLAY THE 
DATA 

15:20 REM 

1540 PRINT " CHOOSE 6 POINTS" 

1550 PRINT 

1550 PRINT "THE CALIBRATION WILL BE STORES IN 'CAL.DAT'" 
1560 REM 

1570 REM Begin sampling loop 

1580 REM 

1600 FOR J < 1 TO 1 

1610 PRINT 

1630 FOR JJ = 1 TO 6 

1631 INPUT "INPUT THE CALIBRATION PRESSURE"; CAL(JJ) 
1632 INPUT "PRESS “ENTER” TO START MEASUREMENT"; MOVES 
1640 FOR II = 1 TO 10 

1650 ROUT = 1 

1660 CALL OUTPUT (RELAY .MUX.01, ROUT) 

1670 CALL MEASURE (DMM.01, VOLTS) 

1520 PA(IL, JJ) = VOLTS 

HEG NEXT 11 

1700 IF JJ = 6 THEN 1740 

moO NEAT JJ 

1740 REM 

1750 REM DISPLAY THE SAMPLE DATA 

1760 REM 

1780 POR IS= 1 TO 10 
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1 


7 


9 O P R I N T U S I N G 


FORMATS ;IS,PA(IS,1),PA(IS,2),PA(IS,3),PA(IS,4) ,PA(IS,5),PA(I 


S,6) 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1920 
1930 
1940 
2000 
2010 
2020 
2055 


PP(J, 


2070 
2080 
2090 
2107 
2102 
2106 
2107 
2108 
2109 
2210 


NEXT IS 

REM 

REM AVERAGE THE DATA 

REM 

FOR JA = 1 TO 6 

TOTAL — 0 

FOR IA = 1 TO 10 

TOTAL = TOTAL + PA(IA, JA) 

NEXT IA 

AVERAGE = TOTAL / 10 

P(JA) = AVERAGE 

NEXT JA 

PRINT 

PRINT "THE AVERAGE ARE: " 

FOR JD = 1 TO 6 

PP(J, JD) = P(JD) 

NEXT JD 

PRINT USING FORMATS; J; PP(J, 1); PP(J, 2); PP(J, 3); 
4); PP(J, 5); PP(J, 6) 

PRINT 

INPUT "DO YOU WANT RE-MEASURE AGAIN ? (Y / N)"; C$ 
IF C$ - "y" THEN 1580 

REM*** STORE DATA BEFORE NEXT SAMPLE*** 

OPEN "A:\CAL.DAT" FOR OUTPUT AS #2 

FOR ID = 1 TO 6 

PRINT #2, USING FORMATS; ID; PP(J, ID); CAL(ID) 

NEXT ID 

CLOSE #2 

NEXT J 
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APPENDIX C 


$storage:2 
Sdebug 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkakk 
k k k k 


x THIS PROGRAM CONVERTS THE VOLTAGE OF TRANSDUCER INTO 
PHYSICAL * 
* PRESSURE, VELOCITY,YAW ANGLE AND PITCH ANGLE. THOSE DATA ARE 

* 
* USED FOR PLOT PROGRAM LATER. 

* 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
k k k k * 


CHARACTER*14 FNAME 
CHARACTER*14 NAME 
CHARACTER*14 OUTFILE 
CHARACTER*2 A(50) 
CHARACTER*80 ST 
REAL K,INTR 
INTEGER COLS,RWS, DTPTS 
BETA A/”017, ”02”, 037,704”, ”05”, ”06”, ”07”, ”08”,”09”, 
^ 10" IR 114” 116517”, : 
7797 57207 727“ 1231 24564127" 
, 7 , 7 
7 , , / 


w * “ ~ 
= 
> 
Ul 
= 


128? E O 27” 254275 Jmm 
137! 7:7 397 ар Aa 143 7447” 
PAG AA “AS 1107775077 
WRITE (*,’(A\)’) ” # OF COLUMNS IN THE GRID 
(LEFT/RIGHT) — ' 
READ (*,*) COLS 
WRITE (*,’(a\)’) 7 4 OF DATA POINTS IN A COLUMN 
(UP/DOWN) = ' | 
READ (*,*) RWS 
WRITE (*,’(A\)’) ’ DATA FILE NAME? (A: FILEXXXX.EXT): ' 
READ (*,’(A14)’) NAME 
WRITE (*,’(A\)’) ’ Initial ambient pressure [in. Hg]: 


+ Æ + H * 


= = = = = mm 
= 
сә 
= 


READ (*,*) PI 

WRITE (*,'(AW)') ' Final ambient pressure [in. Hg]: ' 

READ (*,*) PF 

WRITE (*,’(A\)’) " Initial temperature [deg F]: ' 

READ (*,*) TI 

WRITE (*,’(A\)’) ' Final temperature [deg F]: ' 

READ (*,*) TF 

WRITE (*%,” (AN) ”) ’ Tunnel calibration factor, K: ' 

READ (*,*) K 

WRITE (”,”(AV)”) ” Slope from pressure calibration 
curve (cm/v) ” 
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READ”) EQEE | 
WRITE (*,’(A\)”) * Intercept from pressure calibration: 


READ (*,*) INTR 
WRITE (*,’(A\)’) * Tunnel delta-p for test-section q: 


READ (*,*) OMIFAC 
* COMPUTE THE AVERAGE ATMOSPHERIC PRESSURE GIVEN inHg CONVERT 


TO psf 
PATM=(PI+PF) *35.3631 
R=1716.5 

C ****** "E" is blockage correction *****XXk*kXkXXXX 
Е-0, 0123 


C ****** COMPUTE TEST SECTION AVERAGE TEMPERATURE 
De Ta TR) 72544160 
C ****** TEST SECTION DENSITY 
RO=PATM/ (R*T) 
DTPTS=RWS*COLS 
* OPEN A NEW FILE TO STORE THE REDUCED DATA 
write(*,’(A\)’)’ Output data file? [A:FILENAME.EXT]: ' 
read(*,^(a14)") outfile 
OPEN (2, FILE=outfile,STATUS=’ new’ ) 
WRITE (2,900) DTPTS 
900 FORMAT(I5) 
k OPEN A SEQUENTIAL OF DATA FILE * BE SURE THE FILE IN PPROBE 
HAS 6 ELEMENTS 
DO 20 I=1,COLS 
NAME (9:10)=A(I) 
FNAME=NAME 
OPEN (1, FILE=FNAME) 
READ(1,100,END=20)ST 


100 FORMAT (A65) 
15 . READ(1,1000,END=30)NO,X,Y,V1,V2,V3,V4,V5,BETA 
1000 FORMAT (12,F7.2,F6.2,5F9.3,F8.2) 
* CONVERT THE VOLTAGE TO PRESSURE IN LBF/FT**2 kk 


1cm/H20-2.0461 PSF xxx 
P1=calvp(V1,SLOPE, INTR)*2.0461+PATM 
P2=calvp(V2,SLOPE, INTR)*2.0461+PATM 
P3=calvp(V3,SLOPE, INTR) *2.0461+PATM 
P4=calvp(V4,SLOPE, INTR) *2.0461+PATM 
P5=calvp(V5,SLOPE, INTR) *2.0461+PATM 


IF ((P1-(P2+P3)/2) .LT. 0.0) THEN 


* AN ATTEMP TO ACCOUNT FOR THE DELTA P NOT EQUAL TO 
ZERO * * * k k k k k * k k * * * 


* P1=ABS(P1) 
* P3-ABS(P3) 
* P2<ABS (P2) 
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* CALCULATE THE PITCH ANGLE IN DEGREES 
P-(P4-P5)/(P1-P2) 
ALPHA=FPITCH (P) 
272 LCULATE THE VELOCITY IN FT/SEC 
YSLOP=FYSLOP (ALPHA) 
VELM=SORT ( (2*YSLOP* (ABS ((P1-P3))))/(RO*K)) 
VEL=VELM* (1+E) 
* CALCULATE THE LOCAL DYNAMIC PRESSURE 
ОМ1=ОМ1РЕАС*2.0461/К 
QM=RO*VEL**2/2. 
O1=0M1* (1+2*E) 
O=0OM* (1+2 *E) 
CALCULATE THE YAW ANGLE IN DEGREES 
CG kkkkkkkkkk beta0 is tunnel cal correction ***** 
beta0 = 4. 
YAW-FYAW(BETA - betaO) 


X 


C ***** Calculate the velocity components 
betar = yaw*0.017453 

C *****  alphaO is tunnel cal correction ***** 
alphaO = -2. 
pitch = alpha-alphao 
alphar = pitch*0.017453 
vely = vel*sin(alphar) 
velx = vel*cos(alphar) *sin(betar) 


+ 


CALCULATE THE TOTAL PRESSURE IN LBF/IN**2 
PTC=FPT (ALPHA) 
PT1=P1-Q*PTC 
PT=PT1/144. 
CPT=(PT1-PATM-Q1)/01 


IF (CPT .LT. -3.0) THEN 
CPT=-3.0 
ENDIF 


* CALCULATE THE STATIC PRESSURE IN LBF/IN**2 
Р51=РТ1-О 
PS-PS1/144. 
CPS-(PS1-PATM)/Q1 
GO TO 1190 
ELSE 


* CALCULATE THE PITCH ANGLE IN DEGREES 
P-(P4-P5)/(P1-P2) 
ALPHA=FPITCH (P) 

* CALCULATE THE VELOCITY IN FT/SEC 
YSLOP=FYSLOP (ALPHA) 
VELM=SORT ((2*YSLOP*(((P1-P2))))/(RO*K)) 
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VEL=VELM* (1+E) 
* CALCULATE THE LOCAL DYNAMIC PRESSURE 
QM1=QM1FAC*2.0461/K 
QM-RO*VEL**2/2. 
Q1=QM1*(1+2*E) 
Q=QM* (1+2*E) 
CALCULATE THE YAW ANGLE IN DEGREES 
C *****kk*kkk betaO is tunnel cal correction *****x 
beta0 = 4. 
YAW-FYAW(BETA - betaO) 


% 


C ***** Calculate the velocity components 
betar = yaw*0.017453 

C *****  alphaO is tunnel cal correction *****x 
alphaO Zo: 
pitch alpha-alpha0 
alphar = pitch*0.017453 
vely = vel*sin(alphar) 
velx = vel*cos(alphar) *sin(betar) 


* 


CALCULATE THE TOTAL PRESSURE IN LBF/IN**2 
PTC-FPT(ALPHA) 
PT1=P1-Q*PTC 
PT=PT1/144. 
CPT=(PT1-PATM-Q1)/01 
* CALCULATE THE STATIC PRESSURE IN LBF/IN**2 
PS1=PT1-Q 
PS-PS1/144. 
CPS=(PS1-PATM) /Q1 
ENDIF 
1 1 9 0 
VRITE(2,2000)X,Y,VEL,velx,vely,YAV,pitch,bt,cpt,ps,cps 
2000 FORMAT(11F10.4) 
GO TO 15 
30 CLOSE(1) 
20 CONTINUE 
CLOSE (2) 
STOP 
END 
KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKAKk 
* THIS FUNCTION CONVERTS THE VOLTAGE TO PHYSICAL PRESSURE 
FUNCTION calvp(X,SLOPE, INTR) 


REAL INTR 
calvp=X*SLOPE+INTR 
END 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKAK 
x THIS FUNCTION CALCULATES THE PITCH ANGLE 
FUNCTION FPITCH(X) 
FPITCH=3.759+53.7568*X-1.3085*X*x*2-1.6583*X*x*3 
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* -0.8061*X**4+16.5115*X**5 
END 
KKKKKKAKKAKKAKKKKAKKKKKKKAKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKAKKKKKKKAKAK 
* THIS FUNCTION CALCULATES THE VELOCITY PRESSURE COEFFICENT 
FUNCTION FYSLOP(X) 
IF(X.LT.-10) THEN 
FYSLOP-0.981-0.0102*X-3.000E-4*X**2-2.500E-6*X**3 
ELSE IF((X.GE.-10).AND. (X. LE. 10) ) THEN 
FYSLOP-0.98-0.006*X-2.000E-4*X**2 
ELSE 
FYSLOP=0.9801-0.0035*X-1.143E-4*X**24+5.833E-6*X**3 
END: IF 
END 
kkkkkkkikkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
* THIS FUNCTION CALCULATES THE YAW ANGLE 
FUNCTION FYAW(X) 
C***kkk** NOTE: 180 deg indicates zero yaw 
FYAW=180-x 
END 
kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkik 
* THIS FUNCTION CALCULATES THE TOTAL PRESSURE COEFFICIENT 
FUNCTION FPT(X) 
IF(X.LE.-30) THEN 
FPT--0.01 
ELSE IF((X.GT.-30).AND.(X.LT.-20)) THEN 
FPT=0.02+1.00E-3*X 
ELSE IF((X.GE.-20).AND. (X.LE.30)) THEN 
FPT=0 
ELSE 
FPT=0.03-1.00E-3*X 
END IF 
END 
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